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Abstract. In this paper we want to revisit results of Dahmen and 
Micchelli, |13 ) .|14 j . |15j . which we reinterpret and make more precise. 
We compare these ideas with the work of Brion, Szenes, Vergne and 
others 0, E3, OH, |55). 
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1. Introduction 

The main purpose of this paper is to revisit a celebrated Theorem of Dah- 
men and Micchelli ^H] which we will state and prove again in a somewhat 
stronger form. 

The theorem is on the following settings: 

Start from a finite list X := {a±, . . . ,ajv} of non zero vectors a, G R s . 
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If X spans M s , from X one builds an important function for numerical 
analysis, the box spline implicitly denned by the formula: 

n rX rX N 

(1) / f(x)B x (x)dx := ... fCS2t i a i )dt 1 ...dt N , 

JR S JO JO ~l 

where f(x) is any continuous function. 

If is not in the convex hull of the vectors <Zj then one has a simpler 
function Tx (x) , the multivariate spline cf [EiJ characterized by the formula: 



(2) / f(x)T x (x)dx= [ /(VtiOi)dt, 
Jrs Jr* ~[ 

where f(x) has compact support. 

Both Bx and Tx have a simple geometric interpretation as functions 
computing the volume of certain variable polytopes. 

If furthermore the vectors a, happen to be integral vectors (and is not 
in their convex hull) one has a third function, now on Z n , important for 
combinatorics, the partition function given by: 

(3) V x (v) := #{(ni, . . . ,n N ) | ^n^a; = v, m € N}. 

One of the goals of the theory is to give computable closed formulas for 
all these functions and, at the same time, describe some of their qualitative 
behavior. 

The main result of the theory is that, these three functions can be de- 
scribed in a combinatorial way as a finite sum over local pieces (see formulas 
El and El). In the case of Bx(x) and Tx{x) the local pieces span, together 
with their derivative, a finite dimensional space D{X) of polynomials. In the 
case of Vx(v) they span together with their translates, a finite dimensional 
space V(X) of quasi polynomials (cf. Definition I9.5JI . 

The theorem we are referring to, characterizes: 

• D(X) by differential equations. 

• V(X) by difference equations. 

In particular Dahmen and Micchelli compute the dimension of both spaces. 
This dimension has a simple combinatorial interpretation in terms of X. 
They also decompose V(X) as a direct sum of natural spaces associated 
to certain special points P(X) in a suitable torus. D(X) is then the space 
associated to the identity. 

Although this theorem originates from the theory of the box spline, nev- 
ertheless it has also an interest in commutative algebra and algebraic geom- 
etry, in particular in the theory of hyperplane arrangements and partition 
functions. 
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Here we have been inspired by the results of Orlik-Solomon on cohomol- 
ogy and the results of Brion, Szenes, Vergne on partition functions. 

In fact a lot of work originated from the seminal paper of Khovanskh, 
Pukhlikov interpreting the counting formulas as Riemann-Roch formu- 
las for toric varieties and of Jeffrey-Kirwan, [22] on moment maps. These 
topics are beyond this paper and we refer to Vergne's survey article |38j . 

In the theory of hyperplane arrangements, X is viewed as a set of linear 
equations and the object of study is the arrangement of hyperplanes defined 
by these equations. 

Due to the somewhat high distance between these two fields, people work- 
ing in hyperplane arrangements do not seem to be fully aware of the results 
on the Box-spline. 

On the other hand there are some methods which have been developed in 
this latter theory which we believe shed some light on the space of functions 
used to build the box spline. Therefore we feel that this presentation may 
be useful to make a bridge between the two theories. Thus this paper has 
been organized partly as a research paper and partly as a survey of some 
relevant aspects of these two theories. 

We divide the discussion in three parts. 

In the first on the differentiable case the main new results are determina- 
tion of the graded dimension of the space D(X)( Theorem I6.2|) in terms of 
the combinatorics of bases extracted from X. An algorithmic characteriza- 
tion in terms of differential equations of a natural basis of the top degree part 
of D(X) (Proposition 122(1 from which one obtains explicit local expressions 
for Tx fTheorem l3.6|) . A duality between D(X) and a subspace of the space 
of polar parts relative hyperplane arrangement associated to X (Theorem 
I6.3|) . In this section we also give a simple proof of the Dahmen-Micchelli 
theorem on the dimension on D{X) using elementary commutative algebra 
(Theorem 14.6(1 . 

The second on the discrete case, contains various extensions of the results 
of the first part in the case in which the elements in X lie in a lattice. We 
develop a general approach to linear difference equations and give a method 
of reduction to the differentiable case (Section EJ). We also give an explicit 
formula relating partition functions to multivariate splines ( Theorem 19.22(1 . 
We end with a quick survey of some of the standard applications as for 
example can be found in the book Box splines |18j . 

In the third part we explain the approach via residues (see Theorem lll.lll 
and Formula 160(1 . We finish giving an overview of relations with wonderful 
embeddings. 

Warning: the reader will see that we use Fourier transforms as an es- 
sentially algebraic tool. 
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The fact is that, for the purpose of this theory, Fourier transform is essen- 
tially a duality between polynomials and differential operators with constant 
coefficients. 

As long as one has to perform only algebraic manipulations one can avoid 
difficult issues of convergence and various forms of algebraic or geometric 
(residues) duality is sufficient. 

Thus we usually work with Laplace transforms which avoids cluttering 
the notations with unnecessary i's. 

Our conventions for Laplace transforms are the following. We fix a vector 
space V and set U := V*. We fix a Eucliden structure on V which induces 
Lebesgue measures dv, du on V, U and all their linear subspaces. We set 

Lf(u) := / e-^f(v)dv. 
Jv 

L maps functions on V to functions on U. We have the basic properties, 
when p G U, w G V, writing p, D w for the linear function (p \ v) and the 
directional derivative on V, (and dually on U): 

(4) L(D w f)(u) = wLf(u), L(pf)(u) = -D p Lf(u), 

(5) L(e p f)(u)=Lf(u-p), L(f(v + w))(u) = e w Lf(u). 



Part 1. The differentiable case. 

2. Basic definitions 

It is convenient to take a somewhat intrinsic and base free approach to 
our problems. 

Let us fix an s— dimensional vector space U, let us denote by V its dual and 
fix a list X := {a±, . . . , a^} of non zero elements in V (we allow repetitions in 
the list as this is important for the applications). We identify the symmetric 
algebra S[V] with the ring of polynomial functions on U and sometimes 
denote it by A. 

This algebra can also be viewed as the algebra of polynomial differential 
operators with constant coefficients on V. Similarly S[U] is the ring of poly- 
nomial functions on V, or polynomial differential operators with constant 
coefficients on U. 

Given a vector v G V we denote by D v the corresponding directional 
derivative. This is algebraically characterized, on S[U], as the derivation 
which on an element <fi £ U takes the value {<p\v). 

One can organize all these facts in the algebraic language of Fourier trans- 
form. Let W(V), W(U) denote the two algebras of differential operators with 
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polynomial coefficients on V and U respectively. Notice that, from a purely 
algebraic point of view they are both generated by V U. 

In the first case V is thought of as the space of directional derivatives and 
then we write D v instead of v, and U as the linear functions. In W(U) the 
two roles are exchanged. 

The relevant commutation relations are thus: 

[D v ,<f>] = ((f>\v), [D (f> ,v] = (v\<j>). 
Thus we see that we have a canonical isomorphism of algebras: 

T : W(V) -»• W(U), D v ^-v, 4>^D^. 
One usually writes a instead of J 7 (a) . 

This allows us, given a module M over to consider its Fourier 

transform M as a module over W(U) by a.m := am and conversely. 

2.1. Cocircuits. We come to the first basic definition of combinatorial na- 
ture. The importance of this notion will be clear once we start to study the 
multivariate spline. We assume that X spans V. 

Definition 2.2. We say that a sublist Y C X is a cocircuit, if the elements 
in X — Y do not span V. 

The minimal cocircuits can thus be obtained as follows: 

Fix a hyperplane H C V spanned by elements in X and consider Y := 

{x € X | x ^ H}, it is immediately verified that this is a cocircuit and every 

cocircuit contains one of this type. 

Sometimes we shall express the fact, that a sublist Z C X consists of all 
the vectors in X lying in a given subspace, by saying that Z is complete. 

Thus a minimal cocircuit is obtained by removing from X a complete set 
spanning a hyperplane. 

The set of all cocircuits will be denoted by 8(X). 

2.3. No broken circuits. The second basic combinatorial notion has been 
used extensively in the theory of hyperplane arrangements (cf . , jlU] , , |42| ) . 

Let c := a^, . . . ,ai k 6 X, i\ < ii - ■ ■ < ik, be a sublist of linearly inde- 
pendent elements. 

Definition 2.4. We say that a, breaks c if there is an index 1 < e < k such 
that: 

• i < i e . 

• dj is linearly dependent on a^, . . . , ai k . 

In particular, given any basis b := a^, . . . , Oj s extracted from X, we set: 
B{b) := {a £ X \ a breaks 6} and n(b) = \B(b) \ the cardinality of B(b). 

Definition 2.5. We say that b is no broken if B(b) = b or n(b) = s. 
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Let us denote by B(X) the set of all bases extracted from X. We shall 
consider the map b i— > n(b) ciS cl statistic on B(X). 

2.6. The box spline. Let us recall some points which are standard using 
the form presented in [2] or |2L)| . 

First let us recall some basic facts on splines (cf. |18j ). 

Let C(X) := {^2 aeX ta a I < t a , Va} be the cone of linear combinations 
of vectors in X with positive coefficients. 

We will assume that is not in the convex hull of the vectors in X, i.e. 
that C(X) does not contain lines. 

We have already defined, in the introduction, the two basic functions on 
V, B x (formula and Tx, (formula ©). 

It is best to think of both Tx and Bx as tempered distributions (cf. [39 ). 
Then the definition is valid also if X does not span V. 
Bx is supported in the box 

N 

B(X) := ^2 t i a ^ < < 1, the shadow of the cube [0, 1]^, 
i=l 

generated by X, Tx is supported in C(X). 

Basic example Let X = {a±, . . . , a s } be a basis, d := | det(ai, . . . , a s )\: 
B(X) is the parallelepiped with edges the a{, C(X) is the positive quad- 
rant generated by X. 

(6) B x = d~ l XB{x), T x = d~ l xc{x) 

where, for any given set A, we denote by xa its characteristic function. 

If X = {a±, . . . , cifc}, k < s is only a linearly independent set, we have to 
consider Tx and Bx as meausers on the subspace spanned by X. 

Tx and Bx are functions as soon as X spans V, i.e. when the support pf 
the distribution has maximal dimension. 

These functions have a nice geometric interpretation. 

Let F : R N -> V be defined by F(t u . . . , t N ) := k<H- Then B x (w) 
is the volume of the polytope F~ 1 (w) n [0, 1]^ while Tx(w) is the volume 
of the polytope F^ 1 (w) n [0, oo]^ (with a suitable normalization constant). 

It is useful to generalize these notions, introducing a parametric version 
called E— splines. 

Fix parameters \i := {fi\, . . . ,hn} and define the functions (or tempered 
distributions) on V by the implicit formulas: 

r r 1 f 1 n N 

(7) f(x)B XilM {x)dx := ... e-^ t ^f(Y j t i a l )dt l ...dt N . 
Jv Jo Jo ~1 
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(8) / f(x)T x , u (x)dx = / e^«''"«/(Vt i a i )(li 

Jv - 7r? 

Also these functions have a nice geometric interpretation. 

They represent the integral of e _ ^<=i * <w on the polytope F _1 (w)n[0, 1]^ 
or F -1 (w;) (~l [0, oo]^ (with the same normalization constant). Of course for 
fi = we recover the previous definitions. 

An easy computation gives their Laplace transforms: 



[e-^ y) B x ,^x)dy= f\.. C e~^i ^(M+w)^ . . . dt 
Jv '- Jo Jo 

-n 



N 

I _ g— a— /J, a 



... a + n a 

adX 



(9) 



and 

P f'OO f'OG 

(10) e- {x ' y) T x ,„(x)dy = / .../ e -E i =i*i(^ I ai>+w) dtl . . . d t N 

Jv '- Jo Jo 

= n — 

We have written shortly a := (x, a), for the linear function on U. 

The use of Laplace rather than Fourier transforms is justified by the 
following discussion. 

Define the dual cone C(X) of C{X). 

C[X) := {u € U \(u | v) > 0, Vt)G C(X)}. 

This cone consists thus of the linear forms that are non negative on C(X). 
Its interior in not empty since C(X) contains no lines. 

Proposition 2.7. If T is a tempered distribution supported in C(X) its 
Fourier transform is an analytic function, of a complex variable z = x + iy, 
x G C(X),y G U, on the open set where the real part x lies in the interior 

c{xy> ofc[x). 

Proof. This depends on the fact that, if u G C(X)° we have that e~^"l^ has 
exponential decay on C(X). □ 
In fact for the Tx and all the distributions that we shall encounter we 
will have that their Fourier transform are not only defined in the region 

where the real part x lies in C(X)°, but in fact they extend to meromorphic 
functions with poles on the hyperplanes a% = (or sometimes translates of 
these hyperplanes). 

In the course of this paper we give an idea of the general algebraic calculus 
involving these distributions. Under Laplace transform one can reinterpret 
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the calculus in terms of the structure of certain algebras of rational functions 
(or exponentials) as D— modules. 

Given a £ V, fi € C let us introduce the notation, which will be discussed 
more deeply in Part 2: 

V£/(x) :=f( x )- e -»f(x-a). 

From the expressions of the Laplace transforms one gets that, the box 
spline can be obtained from the multivariate spline by a simple combinatorial 
formula. 

Proposition 2.8. For every subset S C X we set as ■= EaGS a ' an< ^ 
MS : = Eaes then: 

(11) BxA*) = II VpTxJx) = (-l)We-» s TxAx - a s ). 

aGX sex 

Proof. It follows from the basic rule © which gives the commutation rela- 
tion between the Laplace transform and translations. □ 

2.9. The space D(X). The following definition is of central importance in 
the work of Dahmen and Micchelli and in this paper 

Definition 2.10. We define the space D(X) by the condition: 

(12) D(X) := {p | D Y p = 0, VY € £(X), the cocircuits}. 

In this definition we shall assume that p is a polynomial. In fact, due to 
the property that the ideal generated by the elements Dy contains all large 
enough products of derivatives (see 14.4)1 . one can easily see by induction that 
any distribution p satisfying 1)12)1 . is already a polynomial (cf. |26j). 

We will see later a generalization with parameters fi of these equations and 
also a discrete analogue, using difference operators rather than derivatives. 

Let m be the minimum number of elements in a cocircuit in X, assume 
m > 2, m is characterized by the property that the basic space D(X) con- 
tains all polynomials of degree < m. We shall see in Part 4, that m controls 
also the smoothness of the splines that we have introduced. 

2.11. Two basic modules in correspondence. The theory of the Laplace 
transform tells us how to transform some basic manipulations on distribu- 
tions as an algebraic calculus. In our setting this is best seen introducing 
the following two D— modules in Fourier duality: 

The first is the D— module T>x '■= W(V)Tx generated, in the space of 
tempered distributions, by Tx under the action of the algebra W(V) of 
differential operators on V with polynomial coefficients (for basic facts cf. 

Eg)- 

The second I?— module is the algebra Rx := S[V] [Y\ ae x a_1 ] obtained 
from the polynomials on U by inverting the element dx '■= Yla^x a - This 
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is a module under W(U) and it is the coordinate ring of the open set Ax 
complement of the union of the hyperplanes of U of equations a = 0, a € X. 

Theorem 2.12. Under Laplace transform, T>x is mapped isomorphically 
onto Rx ■ In other words we get a canonical isomorphism of T>x with Rx 
as W(U)-modules. 

T>x is the space of tempered distributions which are linear combinations 
of polynomial functions on the cones C(A), A C X a linearly independent 
subset and their distributional derivatives. 

Proof. The injectivity of the Laplace transform on T>x is a standard fact 

To see the surjectivity, notice that under the Laplace transform, by def- 
inition and formula QlOj). the image of T>x under the Laplace transform is 
the smallest -D-module containing d x . Since Rx contains d^ 1 it suffices to 
see that d x l generates Rx as a -D-module. To see this first notice that if we 
take linearly dependent vectors ao = Y2i=i a i a i> m X, we can write: 

1 

So, repeating this algorithm, we see that we can write a fraction l/(n? a f l ) 
as a linear combination of fractions where the elements appearing in the 
denominator are linearly independent. In particular Rx is spanned by the 
functions of the form b^ i+1 ), with / a polynomial, B = . . . bt} C 

X a linearly independent subset and fcj > for each i = 1, . . . , t. 

Now it is clear that //(]X k 1 "*" 1 ) lies in the D-module generated by 
1/01^) and 

1 = UaeX-B a 

Ui b i d x 

So, our first claim follows. 

The second follows from our previous discussion once we remark that 
if B = 6t} C X is linearly independent, 1/(1^6$) is the Laplace 

transform of the measure on V which is the push forward of the measure on 
the subspace spanned by the vectors in B given by Xc{B)h L B where Xc(B) is 
the characteristic function of the cone C(B) and hb is the Lebesque measure 
normalized so that the parallelepiped with edges the bi has volume 1. □ 

Let us now introduce a filtration in Rx by D— submodules which we will 
call the filtration by polar order. 

This is defined algebraically as follows. One puts in filtration degree < k 
all the fractions fY\aex a ~ ha , h a > for which the set of vectors a, with 
h a > 0, spans a space of dimension < k. Denote this part of the filtration 
by Rx,k ■ Notice that by our proof of Theorem I2.12( Rx,s = Rx- 



«o 



in; 



E 



Q:; 



o Hj&.j=i 



a; 
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From the proof of Theorem 12.121 one deduces that the corresponding (un- 
der inverse Laplace transform) filtration on T>x can be described geometri- 
cally as follows. 

We cover C(X) by the positive cones C(A) spanned by linearly indepen- 
dent subsets of A C X. We define C(X)k to be the dimensional skeleton 
of the induced stratification, a union of k— dimensional cones. 

Proposition 2.13. T>x,k consists of the tempered distributions inT>x whose 
support is contained in C(X)&. 

2.14. Polar parts. Proposition ^. 13l implies that, the space T>x,s-i is formed 
by all the distributions in T>x which vanish once computed on test functions 
with support in the set of regular points C(X)° := C(X) — C(X) s _i. 

In other words we may identify T>x /T^x,s-i with a space of distributions 
on the open set C(X)°. 

By duality we then see that the top subquotient, Rx / Rx,s-\, has a special 
importance in the theory [3]. Let us define it in a much more general setting 
which will be useful later. 

Let M be an analytic manifold. Z = Z\ U Z<i U • • • U Z m be a divisor union 
of the smooth irreducible components Zi, i = 1, ... ,m. p an isolated point 
in Z\ n Z<i n • • • n Z m with the property that there exist local coordinates 
around p so that the Z\ are locally given each by a linear equation a, = 0. 
A subset of the set of divisors Z± is said to be non spanning in p, if p is not 
an isolated point of their intersection. We then define: 

Definition 2.15. Define the module of polar parts in p, with respect to Z, 
and denote it by Pz, p (or Pz if p is clear from the context) as the quotient 
of the space of germs in p of meromorphic functions with poles in Z modulo 
those function whose polar part is supported in a non spanning subset. 

Pz,p is clearly a module over the ring Z p of germs of differential operators 
with holomorphic coefficients around p. In the local coordinates x\, . . . ,x s 
around p, where Zi has an equation aj = 0, i linear this space can then 
be identified with the module Px '■= Rx / Rx,s-i with X = {a±, . . . ,a m }. 
Notice the under this identification the operators induced from Z p coincide 
with the algebra of operators W(U) where U = T p X. 

2.16. Basic modules. All the modules over Weyl algebras which appear 
are built out of some basic irreducible modules, in the sense that they have 
finite composition series in which these basic modules appear. It is thus 
useful to quickly discuss these basic modules. Let us take W(V), the differ- 
ential operators on V. The most basic module on W(V) is the polynomial 
ring S[U}. It is the cyclic module generated by 1 and the annihilator ideal 
of 1 is generated by all the derivatives D v , v G V. Its Fourier transform can 
be identified with the module of distributions supported in € U, a cyclic 
module generated by the Dirac distribution 5q. 



12 



C. DE CONCINI AND C. PROCESI 



Given any point p £ U, consider the 1-dimensional S[V] module C p given 
by evaluation at p and the induced W(U) module N p := W(U) ®s[V] ^p- 
N p is easily seen to be irreducible and free of rank 1 as module. 

In fact, in the language of distributions N p is identified to the PF(T^) 
submodule generated by the Dirac distribution 5 p . 

We shall need the following (easy) and standard fact ( 

Lemma 2.17. 1) Given a W(U) — module M and a nonzero element u € M, 
if f v = f{p) v ^f £ SIX], then u generates a submodule isomorphic to N p . 

2) Given linearly independent vectors Ui satisfying the previous hypothe- 
ses, the submodules that they generate form a direct sum. 

3) The module N p has as characteristic variety the cotangent space at p. 

The only use we make of part 3) is its consequence that, for distinct points 
p, q, the corresponding modules N p , N q are not isomorphic. 

Also for a linear subspace A C V or a translate A + v one can define the 
irreducible module Na generated by the Dirac distribution 5a given by: 



The annihilator of 5a is generated by the elements u G A 1 - C U vanishing 
on A (for A + v we have the elements u — u(v), u G A- 1 -) and the elements 
D x , x £ A. the fact that Na is irreducible and the previous elements gener- 
ate the annihilator of 5a can either be verified directly or by remarking that 
Na can be obtained by twisting the polynomial ring by an automorphism 
(a partial Fourier transform) of W(V) defined as follows. One decomposes 
V = A®Z for some complement Z, then one has that W(V) = W(A)<8>W(Z) 
and apply Fourier transform to the factor W{Z). 

In explicit coordinates X{ where A := {xi = 0, i < k} we have that 
the partial Fourier transform is the identity on the variables Xj and their 
derivatives when j > k, while 



As the polynomial ring is a free rank 1 module over the polynomials we see 
that Na is thus a free rank 1 module over the partial Fourier transform of the 

d 

polynomials, that is the polynomial algebra Pa in the variables — — , Vi < 

OXi 

k, Xi Mi > k. 

The module appears naturally in our theory as follows. Take a ba- 
sis c := {ai, . . . , of A and consider the open cone C(c) generated by c. 
Consider next the distribution 5c(c) given by integration on this cone under 
a translation invariant Lebesgue measure for which the parallelepiped gen- 
erated by {ai, ... , a&} has volume 1. Then one easily sees that the Laplace 
transform of 5c(c) * s Wi=i a J l an d Na appears as the module of polar parts 
for the polynomial ring with the a, inverted. 
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In particular fix some set of linear coordinates ai,...,a s for U. If for 
every subset S C {1, . . . , s} we set As to be the subspace where the variables 
cii = 0, G S we see, with the previous notations, that the ring of Laurent 
polynomials decomposes direct sum 

C[af\ af 1 ] = ®sc{i,...,b}Pa s (— ^ )• 

Notice that Paq{tt-^ — ) has as basis, the monomials FT- Qj 4 with fcj < 

Hies a » 

0, Vi G 5 and > 0, Mi £ S. 

This decomposition, together with the description of As as a polynomial 
algebra, gives the explicit partial fraction decomposition for Laurent poly- 
nomials, which is in any case elementary. 

Following the same ideas one can develop the general partial fraction 
decomposition in the case in which we invert any set of linear equations as 
we shall see in Proposition 13.21 and Remark 13.41 



3. The function T x 

3.1. Local expansion . Let us extend the previous ideas to the parametric 
case. We take a list fi := {fi a \a G X} of complex numbers. We introduce 
the ring Rx,^ '■= S[V][Y\ aeX (a + /x a ) -1 ]. It is clear that we can introduce 
on this algebra a filtration completely analogous to that of R x = Rx,o ■ 

To a basis 6 := {hi, . . . , b s }, from X, we associate the unique point Pb&U 
such that bi(j>b) = —fJ-bi for each i = 1, . . . , s. The set P(X, fi) consisting of 
the points pb as 6 varies among the bases extracted from X is called the set 
of points of the arrangement. 

For generic fi all these points are distinct, while for fi = they all co- 
incide with 0. In the other cases we may have various ways in which the 
points Pb will coincide, depending on the compatibility relations among the 
parameters \i. 

Given p G P(X, fi) we define the subset 

X p := {a G X | a(p) + fi a = 0} (X p = U Pb _ =p {b}.) 

of the elements in X such that the affine hyperplane H a of equation a{p) + 
Ha = contains p. It is clear by definition that, if we restrict to the sub- 
set X p , the points of this restricted arrangement reduce to p. Moreover a 
change of variables a' := a + [i a corresponding to a translation, centers the 
arrangement in 0. 

The divisor Z := U ae xHa satisfies the hypotheses of section 0.141 thus 
we can construct, for each point p G P(X,fi) the corresponding module of 
polar parts, which we denote by Pzp- In order to understand these modules 
we need a preliminary construction. 

The following Proposition allows us to reduce the computation of Tx,u to 
that of the various Tx p , for p G P(X, fi). Let us denote, for a no broken basis 
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b by Ub the class of the element naeb( a ~'~ 1 m ^ ne quotient Rx / Rx,s-i- 
Since, given ao £ b we clearly have that (ao + /x ao ) Ilaefe^ + f^a) 1 £ % jS -i 
we verify easily that: 

fu b _ = f( P )u b VfeS[V}. 

It follows that, unless U{, = (which is not the case), it& generates a sub- 
module isomorphic to _/V p . In fact one has: 

Proposition 3.2. 1) The mapping Rx / Rx,s-i ®p^P(x,n)Pz,p is an iso- 
morphism of W(U) modules. 

2) Each Pz,p is an isotypic component. 

3) Each Pz,p is the direct sum of as many copies of N p as the no broken 
bases b in X p , each generated by a corresponding element Ub- 

4) Tx,n = S P gp(x n) c p eP Tx p , with c p some explicitly computable con- 
stants. 

Proof. We can use a slightly more precise reduction. Assume we have lin- 
early dependent vectors ao = X)i=i a « a «> i n X. If v := fj, ao — Y2i=i a «Mai 
we write: 

1 = -l QQ + Map ~ Ya=1 a i( a i + MaJ 

niu ( a i + ) nf=o (°* + ) 

and then, develop into a sum of k + 1 terms in each of which one of the ele- 
ments ai + /j, ai has disappeared. This allows us to separate the denominators 
with respect to the points in P(X,fi). 

If v = we can write 

Ui=o( a i + M»i) (a + Mao) Ui=o( a i + Ma*) 
We expand and then simplify the numerators with the denominators 
and obtain that the element rii=o( a * + Mai) -1 can be expanded as a lin- 
ear combination of elements of type Ila e6( a « + ^a^)~ h \hi > with b a 
basis. Only the elements with all the hi > give non zero classes in the 
module of polar parts, moreover when hi > from the rules of derivatives 
Ila e&( Qi + Ma*) - ' 1 * i s obtained by applying a monomial in the derivatives 
to na 4 ef>( a i + Mai) -1 - This shows that the map is onto. The modules Pz, p 
are contained in Rx / Rx,s-i and belong to distinct isotypic components by 
the previous lemma on the characteristic variety, hence the direct sum of 
Rx I Rx,s-i- 

A similar discussion allows us to reduce the elements appearing in the 
denominator to no broken sets. Thus the final point to verify, applying part 
2) of Lemma \l. 171 is that the elements u& are linearly independent in Pz, P - 
This we shall show as a consequence of the theory of residues (Theorem 
I11.7JI . completing the proof of 1) and 2). 

By (flO|) we see that Ilaex ( a ~ a (p))~~ 1 * s t ne Laplace transform of e p Tx p - 
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Thus we are reduced to showing that, for suitable constants c p we have: 

( is ) n ^7 = £ ° p n = e -^m- 

aex ^ a P eP(X,|t) «ex P ^ a P eP(x,^) aex p ^ 

This follows by induction applying the basic algorithm of separation of de- 
nominators. □ 

Remark 3.3. Notice that the summands Pz,p of Rx/Rx,s-l are the spaces 
of generalized common eigenvalues for the commuting operators induced 
from V (a vector v £V has eigenvalue v(p)). Thus any S[V] submodule in 
Rx I Rx,s-i decomposes canonically into the direct sum of its intersections 
with the various Pz,p- 

There is a similar description of all the pieces Rx k/Rx k— l based on the 
spaces of the arrangement generated by the hyperplanes a = 0, a € X of 
codimension k. 

Proposition 3.4. Rxk/ Rxk-i i- s a direct sum of copies of Fourier trans- 
forms of the modules N\y, as W runs over the subspaces of the arrangement, 
and for given W the sum is made of terms indexed by the no broken bases 
in cC X Pi W 1 - C W 1 - and generated by the class of Ylaec a_1 - 

Proof. Consider, for each W, the map Rxnw- 1 - k c Rx,k —* Rx,k/Rx,k-i- 
By induction the characteristic variety of all irreducible factors in Rxk—i 
and RxnW x fc-l are union of conormal spaces to subspaces A of dimension 
< k - 1. Since Rxnw 1 - kl Rxnw 1 fc— l ^ s a direct sum of irreducibles with 
characteristic variety the conormal space to W, it follows that this map 
induces an inclusion Rxnw 1 - k/Rxnw 1 - fe-l c Rx,k/Rx,k-i- It is easy to 
verify that this map is an isomorphism with the isotypic component of type 
Nw in Rx,k/Rx,k-li reducing to the previous case. □ 

3.5. Local expression for Tx- Formula l)13j) implies immediately that: 

(14) Tx, B (x)= c v^x v {x) 

peP(x, li ) 

This together with formula gives for the box spline: 

BxAx)= Y, c p E (-l) lSl e-^ +p B Xp (x-a s ). 
peP(x,fi) scx-x p 

Of course this is also a reformulation of the identity of Laplace transforms: 

1 _ p-a-^a . . 1 _ --a—Ha 

aex ra peP(x,jj) aex-x p aex p H ' a 

As a consequence of these formulas the essential problem is the determi- 
nation of Tx in the non parametric case. We are indeed ready to state and 
prove the main formula one can effectively use for computing the function 
Tx- 
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We first need some geometry of the cone C(X). Denote by C(b) the 
positive quadrant spanned by the no broken circuit basis b. The positive 
cones C(b) cover the cone C{X) and induce a decomposition of C{X) into 
polyhedral cones which happens to be independent of the order chosen and 
thus of the no broken circuits (this is proved in [U). The points of C(X), 
outside the boundaries of these cones are called regular and they are a union 
of open convex cones called the big cells. 

The remaining singular points are the union of cones C (A) generated by 
subsets A C X which do not span V. We denote by MB the set of no broken 
bases extracted from X. 

Theorem 3.6. Given a point x in the closure of a big cell c we have 
(15) T x (x)= |det(6)|-V, x (-z). 

b | cCC(6) 

where for each no broken basis b, Pb,x(x) is a uniquely defined homogeneous 
polynomial of degree \X\ — s lying in D(X). 

Proof. By the continuity of T x in C(X), it is sufficient to show our claim in 
the interior on each big cell. Thus, by formula © and our discussion above 
we need to show that the identity ((To]) holds, with Pb,x(x) G D(X), in the 
D-module Vx/D x , s -i- 

If we work in the space of polar parts Px = Rx / Rx,s-i we have, by our 
previous discussion that there exist uniquely defined polynomials Pb,x(x) 
homogeneous of degree |X| — s with 

< 16 > 

A beMB - 

In fact this identity already holds in Rx- 

Applying the inverse of the Laplace transform (cf. |H1 ll()j) we obtain iden- 
tity (fTKll . 

To finish notice that Y is a cocircuit we have Wy^YV^x 1 e Rx,s-i- hi 
other words, the polynomials p^x i n formula 1)16(1 lie in D(X). 

□ 

From now on, unless there is a possibility of confusion, we shall write pt, 
instead of pb t x ■ 

Remark 3.7. The polynomials Pb(x), with b a no broken basis, are character- 
ized among the polynomials in D(X), by the further differential equations 
which will be described at a later section (see ()22(0 . In particular we will 
show that the polynomials Pb{x) are linearly independent. 

4. The first theorem 

4.1. The basic equations. In section 2.4 we have introduced the elements 
Dy := Ilxey D x € A, Y C X, thought of as differential operators on V. 
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We have started to see that the space 

(17) D(X) := {p G S[U] | D Y p = 0, VFe £(X)}, 

of polynomials p on 1^, which are solutions of the differential equations, 
DyP = 0, for all the cocircuits Y, plays a fundamental role in the determi- 
nation of T x . In the authors prove that: 

Theorem 4.2. Let D(X) C S[U] be the space of polynomials on V , solutions 
of the differential equations Dyf = 0, as Y runs over the cocircuits in X. 

D(X) is finite dimensional, of dimension the number d(X) of linear bases 
ofV which one can extract from X. 

Since the elements Dy are homogeneous, the space D{X) is also graded 
and it is interesting to compute its dimension in each degree D^X). 

Moreover from general facts we shall see that D(X) is generated, under 
taking derivatives, by the homogeneous elements of top degree N — s. 

As usual one can arrange these dimensions in a generating function, 
Hx(q) '■= dim(Z)fc(X))g fc . We will show (Theorem I6.2|) . by exhibiting 
an explicit basis, that this polynomial is given by the statistic introduced in 
§2.1, i.e.: 

(18) H x (q):= £ q N ^\ 

beB(X) 

In particular the top degree polynomials in D(X) are of degree N — s. They 
appear in the formula of the multivariate spline and thus have a very inter- 
esting geometric interpretation, since they compute the volume of certain 
polytopes associated to X ( 20 ). 

4.3. A remarkable family. In order to obtain the proof of 14.21 and Q18|) . 

we start from the study of a purely algebraic geometric object. 

For notational simplicity let us denote by A := S[V]. 

We want to describe the scheme defined by the ideal Ix of A generated 
by the elements My := n^eY x as ^ rims ° ver all the cocircuits. 

Thus we are interested in the algebra: 

(19) A X :=A/I X . 

We shall soon see that A x is finite dimensional. Formally A x is the coordi- 
nate ring of the corresponding scheme. 

The use of the word scheme may seem a bit fancy. What we really want 
to stress by this word is that we have a finite dimensional algebra (quotient 
of polynomials) which as functions vanish exactly on some point p but at 
least infinitesimally they are not constant. This appears clearly in the dual 
picture which produces solutions of differential equations. 

To warm us up in the proof let us verify that this scheme is supported at 
(which implies that A x is finite dimensional). 
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For this, remark that the variety of zeros of a set of equations, each one 
of which is itself a product of equations, is the union of the subvarieties 
defined by selecting an equation out of each product. Thus what we need is 
the following: 

Lemma 4.4. Take one element Xi from each cocircuit Yi, then the resulting 
set of elements Xi span V (hence they define {0} as subvariety). 

Proof. If, by contradiction, these elements do not span, their complement 
is a cocircuit. Since we selected an element from each cocircuit this is not 
possible. 

□ 

In fact it is convenient to extend the notion as follows. If X = {a\, . . . , a^} 
and fj, := {/xi, . . . , /ijv} are parameters, we can define in general the ideal 
Ix{p) given by the equations Y\ a eY ( a j ~ A*i)> ^ ^ £(^0- This can either 
be viewed as an ideal in S[V] depending on the parameters [i or as an ideal 
in the polynomial ring 5[V][/xi, . . . , /x/v]- 

4.5. The first main theorem. Theorem (|4.2j) follows easily from: 

Theorem 4.6. For all \x the ring Ax(n) ■= A/Ix{^) has dimension equal 
to the number d{X) of bases which can be extracted from X. 

Proof. We will use a standard procedure in commutative algebra. First we 
show that dim(^4x(^)) > d(X) for generic \i. Then, for the special point 
/i = we show that dim(ylx) < d{X). Since the dimension is semicontinuous 
this implies the statement (we will comment after about the meaning of this 
theorem in terms of flatness and the Cohen Macaulay property). 

As a first step we claim that, for generic ii, the ideal Ix(^) defines d(X) 
distinct points, (where d(X) is the number of bases extracted from X). 

This implies of course that dim(Ax(^)) > d(X) where equality means 
that all the points are reduced. 

Set theoretically, the variety described by all the equations given is the 
union of the varieties described by selecting, for every Y € £{X) a cocircuit, 
an element a € Y and setting the equation a — fi a = 0. 

From Lemma 1)4.3(1 , the resulting list of vectors , . . . , aj M (extracted 
from all the cocircuits) generates the space, thus the equations — fj,^ = 
are either incompatible or define a point of coordinates Hi- in some basis 
extracted from X. 

Conversely, given such a basis every cocircuit must contain at least one 
of the elements hence the point of equations = fj-ij must be in the 
variety. 

Now it is clear that, if the /Xj are generic (i.e. do not satisfy some linear 
compatibility equations) these d{X) points are all distinct. Proving the first 
inequality. 
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We now prove that dim(Ax) < d(X). To do this we proceed by induction 
on the cardinality of X and dim(V). 

Take y G X and set Z := X - {y}. Set V := V/Cy and denote by 
7r : V — > V the quotient homomorphism. 

Also, for any sublist B C Z set B equal to the list of non zero vectors in 
V/Cy which are images of vectors in B. 

Define finally A := A/(y) = S[V/Cy] = S[V}. 

(1) The image of Ix in A is the ideal 1-%: 

The image of Ix is generated by the products Ilxey n ( x ) wnere Y £ £P0 
and y $.Y. Thus X — Y generates a proper subspace of V containing y. This 
clearly implies that ir((X - Y)) C V and that ~Z n n((X — Y}) = X - Y. 

Therefore the image ir(Y) = Y is a cocircuit for Z and 

n 7r ( x ) = n w - 

This proves that the image of Ix is contained in 1-% (the other generators 
map to 0). 

On the other hand we know that I-^ is generated by the products Il^gy x 
with (Z-Y) C V/Cy. Setting Y = X — vr _1 ((Z - F)), we immediately see 
that Y G £{X) and n^eY^ = FLey 71 "^) proving that coincides with 
the image of Ix- 

Therefore by induction dim(A/(I x + (y))) = din^A^-) = d(Z). 

(2) The surjective map A yA ^> yAx of A modules (multiplication by 
y followed by the projection px : A — > Ax j factors through the projection 
pz ■ A — > A^ followed by a (surjective) map: 

j : A z -> y^x- 

To see this it suffices to show that, if Y G £(Z) then y • rXrgy x ^ Ix- 
We have two cases. If y G (Z-Y), then Y = X - (Z- Y) and Y[ xeY x G 
/a- 

Ify i (Z -Y), then Xn (Z-Y) = Z-Y and X- (Z-Y) =X-(Z-Y) = 
Y U {y} G so that y • Jl^y s = U xe Yu{ y } x ^ J x- 

In particular this proves that dim(A^) > dim(yAx)- 

Let us now use these facts to estimate dim(Ax). 
Consider the exact sequence 

-► yA x -> A x A/{I x ,y) -> 0. 
By what we have seen we have 

dim Ax = dim A/(I x ,y) + dimy^x < dim(A^) + dim(A z ). 
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By our inductive hypothesis, Az has dimension equal to the number of 
bases of V contained in Z i.e. the number of bases in X not containing y 
among its elements. 

Also by induction dim(A^) equals to the number of bases of Vj (x) con- 
tained in Z. 

Given such a basis 6' = {^(a^), . . . , ^(a^-J}, the set b = {y, a^, . . . ai 3 _ 1 } 
is a basis of V contained in X and containing y as an element. Vice versa, 
given a basis b = {y, , . . . fli s _ 1 } as above, then b' = {ir^a^), . . . , Tr(ai s _ 1 )}, 
is a basis of V j (x) contained in Z. 

Thus dim(Ag-) equals the number of bases of V contained in X and con- 
taining y. 

Summarizing dim^z + dim(A^) = d{X). The theorem is proved. 

□ 

Remark that, as a consequence of the proof, which we will use later, we 
also now have that j : Az — > yAx is bijective. 

As announced we want to discuss the meaning of this theorem for the 
variety Vx given by Ix(p) thought of as ideal in 5[V][/ii, . . . , /ijv]- 

This variety is easily seen to be what is called a polygraph. It lies in 
UxC N and can be described as follows. Given a basis b := {a. L1 , . . . , aj s } of 
V, extracted from X, let a 1 , . . . , a s be the associated dual basis in U. Define 
a linear map if, : C N — > U by: 

n 

ib Oil,---, Hn) ■= z^^j a%3 - 

Let Tft be its graph, then Vx = ^b^b- 

Vx comes equipped with a projection map p to C N whose fibers are the 
schemes defined by the ideals 

Remark 4.7. Theorem 14.61 implies that, p is flat of degree d{X) and Ix{p) 
is the full ideal of equations of Vx- Furthermore Vx is Cohen Macaulay. 

The fact is remarkable since it is very difficult for a polygraph to sat- 
isfy these conditions. When it does, this has usually deep combinatorial 
implications (see |28])- 

One should make some remarks about the algebras Ax{p) in general. 

First some basic commutative algebra. Take an ideal I C C[a?i, . . . ,x m ] of 
a polynomial ring. C[x±, . . . ,x m ]/I is finite dimensional if and only if the 
variety of zeroes of / is a finite set of points pi,...,pj~. In this case moreover 
we have a canonical decomposition 

C[xi, . . .,x m ]/I = ef =1 C[xi, . . .,x m ]/I Pi 

where, for each p G {pi, . . . ,Pk}, the ring C[a?i, . . . , x m ]/I p is the local ring 
associated to the point p. 
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Let p have coordinates Xi = /Xj, the local ring C[sci, . . . ,x m }/I p is char- 
acterized, in terms of linear algebra, by the property that the elements Xi 
have generalized eigenvalue /ij. Thus the previous decomposition is just 
the Fitting decomposition, into generalized eigenspaces, for the commuting 
operators X{. 

In the case of the algebra Ax(fi) quotinet of S[V] by the ideal Ix{p) 
generated by the elements, n?/gy(y + %) we see that, if for a point p of the 
resulting (finite) variety we have that y(p) + n y ^ 0, in the local ring of p 
the element y + ^ y is invertible and so it can be dropped from the equations. 
We easily deduce that: 

Proposition 4.8. The local component Ax(fJ>)(p) equals the algebra Ax p (n p ) 
defined by the sublist X, p := {x € X \ x{p) + fi x = 0}. 

Furthermore by a change of variables X p (n) := {x + fj, x , x € X p } we can 
even identify A Xp {[J- p ) = A/r p {I Xp )- 

By t p (Ix p ) we mean the ideal Jv translated at p by the automorphism 
of A sending x i— ► x — x(p),Vx € V (hence x i— > x + fi x , Vx S X p ). 

5. Solutions of differential equations 

5.1. Differential equations with constant coefficients. Let us reinter- 
pret Theorem 14.61 in terms of differential equations. 

Let us consider the polynomial ring S[V], its graded dual is S[U\. 
By duality S[V] is identified with differential operators with constant 
coefficients on V and, using coordinates, the duality pairing can be explicitly 

d d 

described as follows. Given a polynomial p(— — , . . . , — — ) in the derivatives 

ox i ox s 

and one q(xi, . . . , x s ) in the variables the pairing 

(20) <*&---sbM* 1 *•» 

is obtained by applying p to q and then evaluating at 0. 

Of course the algebraic dual of S[V] is a rather enormous space of a rather 
formal nature. It can be expressed best as a space of formal power series by 
associating to an element / € (5[V])* the formal expression: 

k=0 

Where now (/ | x k )/kl is a genuine homogeneous polynomial of degree k on 
V. 

The following facts are easy to see: 

Proposition 5.2. (1) Given a vector v G V, the transpose of multipli- 
cation by v is the directional derivative D v in (S[V])*. 
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(2) Given <f> G U, denote by the automorphism of S[V] induced by 
translation x i— > x — x((p) = x — ((f) \ x), Vx € V, and by its trans- 
pose. We have rlf = e - ^'^ /. 

Proof. 1) follows from the fact that for (j) € U D v ((p) = (<p\v). and the chain 
rule. 

As for 2), 

W I e x ) = (/ 1 r^e*)) = (/ | e*-<* ' ) = e~<+ ' (/ | e*). 

□ 

Observe that, if J is an ideal of S[V] defining a subvariety Z C U, we 
have that J) defines the subvariety Z + 0. 

If we have a quotient S'fV]// by an ideal I we get an injection of % : 
(S[V]/I)* — > (5[V])*. The image is, at least formally, the space of solutions 
of the differential equations given by /. We denote by Sol (I) the space of 
C°° solutions of the system of differential equations given by / 

Assume now that S[V]/I is finite dimensional. Denote by {</>i, . . . , (pk} C 
?7 the finite set of points which are the support of /. Take the decomposition 

s[vyi = ®Us[vyi{<t>i) 

where S[V]/I((f>i) is local and supported in Under these assumptions we 
get 

Theorem 5.3. 1) If S[V]/I is finite dimensional and supported at 0. The 
image of i* lies in S[U] and coincides with Sol(I). 

2) If S\V]/I is finite dimensional and supported at a point <p € U, 

Sol( y I) = e {(t ' lx) Sol{T^I). 

3) For a general finite dimensional S\V]/I, 

Sol(I) = ® k i=1 Sol{I(4>i)) 



Proof. 1) is the duality. As for 2), clearly, for any ideal / and (f> we have 
SoI(t^{I)) = r*^(£of(I)) = e^\ x ^Sol{I). It only remains to explain the 
analytic nature of the statement. This follows, even in the stronger setting 
of tempered distributions by the fact the a tempered distribution which is 
annihilated by high enough derivatives is necesarily a polynomial (see 
Prop. 3.1). □ 

Applying this discussion to our setting, we get that, for generic fx, the 
solutions of the differential equations of the ideal Ix{fJ>) are exactly the space 
with basis the functions where the elements 4> run on the d(X) points 
defined previously, while for fi = we have a remarkable space of polynomial 
solutions. The nature of this space is explained in the next section. 
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In the general case we have a mixture between these two extreme cases. 
The set of solutions is a finite set of points and in each such point <j) G U 
we have a subset of X formed by those elements Xi £ X such that 
(<p\xi) = m. 

Then one can see that this point contributes to the solutions with the 
functions e^p where p are the polynomials associated to X^ at 0. 

6. A REALIZATION OF Ax 

6.1. The graded dimension. Our next task is to prove that, the graded 
dimension of the space D{X) is given by: 

Theorem 6.2. 

Hx( q )= £ q N ~ n ®. 

beB(X) 

In order to do this, we want to realize the ring Ax as the S[V] submodule 
Qx , of the space of polar parts Px generated by the class ux , of the element 
dZ 1 . In doing this we shall exhibit a homogeneous basis made by elements 
Ub as b varies among the bases of X, each of degree N — n(b). 

The spaces Rx and Px are naturally graded (as functions). It is conve- 
nient to shift this degree by ./V so that ux has degree and the generators Ub 
have degree N — s. If b is a no broken basis, these will be just the elements 
Ub introduced in section mi when fi = 0. 

With these gradation the natural map ir : A — > Qx defined by ir(f) = fux 
preserves degrees. 

It is clear that Theorem ()6.2|) follows from the following more precise 
result on Qx, which also describes a graded basis for Ax- 

Theorem 6.3. (1) The annihilator of ux is the ideal Ix generated by 
the elements My = ELey x > as ^ runs over the cocircuits. Thus 
Qx — Ax as graded A-modules. 
(2) Given a basis b := {a^, . . . ,aj s } extracted from X, set 

Ub_:=( ]| a)u x - 

aeX-B(b) 

The elements Ub, as b runs over the bases extracted from X, are a 
basis of Qx ■ 

Proof. To see the first part, let us remark that clearly Ix^x = since, if 
Y is a cocircuit Myux = Y\xex-Y x l nes m Rx,s-i, hence it is in the 
module Px- 

Therefore, from Theorem 14.61 it is enough to see that dimQx > d(X). 

We want to proceed by induction on s and on the cardinality of X. 

If X consists of a basis of V, clearly both Ax and Qx are 1-dimensional 
and the claim is clear. Otherwise we can assume that X := {Z, y} and Z 
still spans V. We need thus to compare several of the objects under analysis 
in the case in which we pass from Z to X. 
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Let us consider the ring A/(y), polynomial functions on the subspace of 
U where y vanishes and denote by Z the set of non zero vectors in the image 
of Z (or X) in A/(y). 

As in Theorem (|4,6|) the set Bx of bases extracted from X can be decom- 
posed into two disjoint sets. 

Bz and the bases containing y. This second set is in 1-1 correspondence 
with the bases of V/Cy contained in Z. 

We need several lemmas. 

First we obviously have an inclusion Rz C Rx and also Rz,k C Rx,k, Vfe. 

Lemma 6.4. Rz >s -i = Rz^Rx,s~i and we get an inclusion of Pz into Px- 

Proof. The two statements are equivalent and follow immediately by order- 
ing Z and adding y as last element. Then all the no broken circuit bases 
for Z are also no broken circuit bases for X so that Pz (as module over the 
differential operators with constant coefficients) is a free direct summand of 
Px- □ 
As a consequence of this Lemma let us consider in Px the map multipli- 
cation by y. We have clearly yux = uz thus we obtain an exact sequence 
of A— modules: 

O^K^Qx^Qz^O 
where K = Qx H Ker(y). 

We need to analyze K and prove that dim(i^) is greater or equal to 
the number d y (X), of bases extracted from X and containing y. Since we 
already know that dim(K) + d{Z) < d(X) and d(X) = d(Z) + d y (X) this 
will prove the claim. 

In order to achieve the inequality dim(i^) > d y {X) we will find inside K 
a direct sum of subspaces whose dimensions add up to d y (X). 

Let us first discuss a special case. Assume that Z spans a subspace V of 
codimension 1 in V, so y is a vector outside this subspace. 

We clearly have inclusions: 

Rz c R x y^Rz,k-i C Rx,k Vfc. 

Also the element ux £ Px is killed by y and A/(y) can be identified to 
S[V'}. 

Lemma 6.5. The multiplication by y^ 1 induces an isomorphism between 
Pz and the kernel of the multiplication by y in Px ■ 

Proof. Since y~ 1 Rz,k-i C Rx,k it is clear that the multiplication by y~ l 
induces a map from Pz = Rz,s-i/ Rz,s-2 to Px = Rx, s /Rx,s-i- It is also 
clear that the image of this map lies in the kernel of the multiplication by 

y- 

To see that it gives an isomorphism to this kernel, order the elements of 
X so that y is the first element. A no broken circuit basis for X is of the 
form {y,c} where {c} is a no broken circuit basis for Z. 
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Now fix a set of coordinates xi,...,x s such that x\ = y and X2, ■ ■ ■ , x s is 
a basis of the span of Z. Denote by di the corresponding partial derivatives. 
We have that in each summand C[di,..., d s ]u y£ the kernel of multiplication 
by x\ coincides with C[e?2, . . . , d s ]u V: c- 

The claim follows easily since C[^2, • • • , d s ]u V) c is the image, under y^ 1 of 
C[d 2 , ■ ■ ■ ,d s ]uc- □ 

Notice that y~ 1 uz = Ux, so that the following lemma is immediate 

Lemma 6.6. Under the previous hypotheses, the multiplication by y^ 1 in- 
duces an isomorphism between Qz and Qx- 

Let us now pass to the general case. Consider the set S y {X) of all complete 
sublists of X which span a subspace of codimension 1 not containing y. 

For each Y € S y {X) we have, by Lemma 16.51 that the multiplication by 
y^ 1 induces an inclusion iy : Py — > Pyu{y} ~^ Px with image in the kernel 
Ker(y) C Px- Thus we get a map 

g := ®YeS y {x)iY ■ ®YeS y (X)Pv -> Ker(y) 
We then claim that: 

Lemma 6.7. g is an isomorphism of ©yg5(x)-Py onto Ker(y). 

Proof. As before order the elements of X so that y is the first element. A 
no broken circuit basis for X is of the form {y,c} where {c} is a no broken 
circuit basis for Y := X n (c). 

By construction Y € S y {X) and we obtain the direct sum decomposition 

Px = ®YeSy(X) P Yu{y}- 

Clearly 

Ker{y) = ®YeS(X)pYu{y} H Ker{y) 
and, by Lemma |6.51 for each Y € S y (X), i y gives an isomorphism of Py 
with Pyvj{ y } H Ker(y). This proves the lemma. 

□ 

In order to finish the proof of Theorem 16.^1 notice that by Lemma ldfil and 
Lemma lfT7| we get an inclusion of ®y^s v (X)Qy into K = Qx H Ker(y), so 
that dimQx > dimQ^ + Syg<Sj / (x) dimPy = d(-X"). 

This gives the required inequality and implies also that we have a canon- 
ical exact sequence of modules: 

(21) -» ®y eSy{X) QY -^Qx^Qz^O. 

It remains to establish the second part of Theorem We prove it by 
induction, ordering X so that y is the last element. Let us write Bx as 
the disjoint union of Bz and the set B^ of bases containing y. Also set 
Ux = {ub | b € Bx} and write it as the dijoint union of U' x = {ub \ b £ £>z} 
and its complement = {u& \ b € Sj^ }. Remark that by our definitions, 
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yli'x = Mz while U' x is the disjoint union Uy g s (x^y^y)- By induction 
and formula (|21|) this clearly implies that Ux is a basis of Qx- D 

It follows from our theorem that in top degree Qx has as a basis the 
elements Ub as b runs over the set of no broken bases. 

In the parametric case, we have the decomposition d^ 1 = c v&~x\ ( c ^- 
(fT3|0 , We set Qx(fJ-) equal to the iS[y]-module generated by ux and for 
each p € P(X,fj,), Qx{p) equal to the S^J-module generated by ux p - 

Proposition 6.8. 1) 

Qx(f±) = ® v &p(x p ^ p )Qx p {}i p )- 

2) QxifJ-) is isomorphic to Ax(^) and the above decomposition coincides 
with the decomposition of Ax(fJ-) into its local components. 

Proof. The first part follows immediately from Remark 13,31 

As for the second since the annihilator of ux contains Ix„ , we have a 
map Ax(^) — * QxilA which by Fitting decomposition induces a map of the 
local summands. On each such summand this map is an isomorphism by 
the previous theorem since we can translate p to 0. From this everything 
follows. □ 

We can easily prove as corollary a theorem by several authors, see P, 

mm- 

Corollary 6.9. Consider in S[V] the subspace V(X) spanned by all the 
products My ■= TlrreY x i Y C X such that X — Y spans V . Then V(X) is 
in duality with D(X). 

Proof. Multiply by d x . V(X)d x 1 is spanned by the polar parts Y\ xeZ x" 1 , 
Z C X such that Z spans V. We have seen that this space of polar parts 
maps isomorphically to its image into Px and its image is clearly Qx- This 
proves the claim. □ 

Remark 6.10. The last theorem we have proved is equivalent to saying that, 
in the algebra Rx the intersection d x l A n Rx, s -i = d^Ix- 

It is of some interest to analyze the deeper intersections d x S[V]C\Rx,k '■= 
Uk- We will sketch this point which uses the structure of the filtration by 
polar degree as explained in |22j . 

If X does not necessarily span V we define Ix as the ideal generated by 
the products Ilxey x where Y C X is any subset such that the span of 
X — Y is strictly contained in the span of X. We set AxiV) = A/Ix- 

Of course if we fix a decomposition V := {X) © T we have Ax(V) = 
A X ®S[T}. 

Consider the set of all subspaces W spanned by elements of X (including 
the space {0} spanned by the empty set). We call any such subspace a 
subspace of the arrangement generated by X. Set L = d^A, X\y := X(l W, 
and Lk = Ln Rx,k- We have 
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Theorem 6.11. For each k we have that Lk/Lk-i is isomorphic to the 
direct sum ®Ax w (V) as W varies on the subspaces of dimension k of the 
arrangement. 

Proof. Under the map A — > L given by / — > d x x f, we know that the submod- 
ule Ix maps isomorphically onto L s -x. So L s _i is spanned by the elements 
dy 1 where Y = X\y for the hyperplanes W of the arrangement. By (a small 
generalization of) Lemma 16.41 we know that the filtration by polar order in 
Rx induces the same filtration in Ry thus by induction we have that the 
graded associated to <iy 1 5[y] is the direct sum of all the pieces Ax w as W 
varies on the subspaces of dimension k of the arrangement generated by Y . 
This follows from Proposition 13.41 □ 

The previous theorem gives an interesting combinatorial identity once we 
compute the Hilbert series oid x l S[V] directly or as sum of the contributions 
given by the previous filtration. 



Remark 6.12. In Fourier transform the previous discussion can be translated 
into an analysis of the distributional derivatives of the multivariate spline, 
i.e. into an analysis of the various discontinuities achieved by the successive 
derivatives on all the strata of the singular set. 



In the applications to the box spline it is interesting, given a set X of 
vectors which we list is some way, to consider for each k > the list X k in 
which every element a € X is repeated k times. Let us make explicit the 
relationship between Hx(q) and H X k(q). 

First the number of bases in the two cases is clearly related by the for- 
mula d(X k ) = d(X)k s , to each basis b := {b\, . . . , b s } extracted from X we 
associate k s bases b(h\, . . . , h s ), indexed by s numbers hi < k expressing the 
position of the corresponding bi in the list of the k repeated terms in which 
it appears. 

Now it is easy to see that: 




n(b(h\, . . . , h s )) = k(n(b) 



s) + hi H \-h s . 



Thus we deduce the explicit formula: 




b&B x hi,...,h 3 
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6.13. More differential equations. Our next task is to fully characterize 
the polynomials Pb(x) appearing in formula (|16j) by differential equations. 
In Theorem 13.61 we have seen that these polynomials lie in D(X). 

For a given no broken circuit basis b, consider the element Db := na0 a - 

Proposition 6.14. The polynomials pb satisfy the system 



(22) DbPcixx, . . . ,x s ) 



1 if b = c 
if b^c 



In particular the polynomials pb for b G MB are linearly independent and 
characterized, in D(X)n-s by the equations \2ty) . 

Proof. This follows from the identity DbUx = the linear independence 
of the elements u& and the fact that the dimension of D(X)x- s is the car- 
dinality of MB. □ 

As a consequence, using formula (|15l) . we can characterize by differential 
equations the multivariate spline Tx{x) on each big cell c as the function T 
in D(X)n-s satisfying the equations: 



(23) D b T 



Idet®!" 1 if ccC(6) 
otherwise 



We have identified D(X) to the dual of Ax and hence of Qx- The basis 
Ub, we found in Qx, defines thus a dual basis u- in D(X). 

Corollary 6.15. When b G MB we have, as polynomials u- = pb 

Proof. This follows from the previous Proposition and the definition of the 
duality given in formula (|20|) . □ 



We shall use a rather general notion of the theory of modules. Recall that 
the socle s(M) of a module M is the sum of all its irreducible submodules. 
Clearly if N C M is a submodule s(N) C s(M) while for a direct sum 
s(Mi M 2 ) = s(Mi) s(M 2 ). If M is finite dimensional s(M) ^ so that, 
for a non-zero submodule N we must have iV n s(M) ^ 0. 

Proposition 6.16. The socle of the S[V]-module Qx coincides with its top 
degree part, with basis the elements Ub- 

Proof. The socle of the algebra of constant coefficients differential operators, 
thought of as a module over the polynomial ring, is clearly generated by 1. 
It follows that the socle of Px (as S'fyj-module), has as basis the elements 
Ub- Since Ub G Qx we have s(Qx) = s(Px), the claim follows. □ 

Theorem 6.17. D(X) is spanned by the polynomial pb, as b runs over the 
set of no broken bases, and all of their derivatives. 
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Proof. Everything follows once we observe: 

1) D(X) is in duality with Ax and so with Qx ■ 

2) The polynomials pb, b € MB are in duality with the elements Ub- 

3) The orthogonal of a proper submodule N of D(X) is a non-zero sub- 
module, thus it must intersect the socle of Qx in a non trivial subspace. In 
particular N cannot contain all the elements pb- □ 

We finish this section observing that we can dualize the sequence (1211) to 
get an exact sequence 

_» D(Z) ^ D{X) — > ® Y eS y (X)D(Y) - 0. 

Here, since the map Qx — * Qz — > is given by multipliation by y, the 
inclusion — > D(Z) — > D(X) is the inclusion of D y (D(Z)). 

Part 2. The discrete case 

7. The discrete case 

7.1. The case X in a lattice. Splines are used to interpolate and approx- 
imate functions. For this purpose it is important to understand the class of 
smoothness of a spline. In the case of the box-spline it is easy to prove (see 
CHI): 

Let m be the minimum number of elements in a cocircuit in X, assume 
m > 2, then Bx is of class C m ~ 2 . 

Thus, provided we choose the list X in a suitable way, we can achieve any 
finite level of smoothness required. 

A particularly useful case is when we have chosen a lattice A C V spanning 
V such that each vector in the list X lies in A. We assume that the Lebesgue 
measure on V is normalized in such a way that a fundamental domain for 
A has volume 1. 

In this case, if C s i ng {X) denotes the set of singular vectors of the cone 
C(X), the set of all translates \J\ e f V C s i ng {X) is called the cut region. It is 
a union of a finite number of s — 1— dimensional bounded poly topes and all 
their translates. 

From formula (|lljl. and the fact that D(X) is stable under translation we 
obtain (see also [T5]): 

Proposition 7.2. The complement of the cut region is a union of all trans- 
lates of a finite number of cells, each an interior of a (compact) polytope. 

Over each such cell the functions Bx(x — A), A £ A are polynomials in 
the space D(X). 

One easily proves the fundamental fact (see [T8]): 
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Theorem 7.3. If X spans V , the translates of Bx form a partition of unity: 



(24) l = Y,Bx(x-X). 

AeA 

Proof. If X is a basis, Bx is the characteristic function of the parallelepiped 
with basis X divided by its volume and it easily follows that the identity 
(|24*|) holds outside a set of measure 0. In general one can use the iterative 
description of the box spline obtained by integration: 

Bx,v(x) = / B x (x -tv)dt, 
Jo 

a formula deduced immediately from the definition (^Q). 



r-l r i N 

Jo 



„1 N „i „ 

/ f(y tjaj)dti . . . dtx = / / f(x + tv)Bx(x)dxdt 
Jo i=1 Jo Jr s 



f(x)B x (x-tv)dxdt= [ f(x)([ B x (x -tv)dt)dx. 

Jr s Jo 

□ 

Given a function /(A) on the lattice A define 

B x *f(x) :=J2 B x(x-X)f(\) 

AeA 

The space of all functions obtained by this procedure is called the cardinal 
spline space and the space of polynomials D(X) is characterized by the 
property of consisting exactly of the polynomials lying in the cardinal spline 
space (see |18|). An extensive portion of the literature on the box spline 
is devoted to understanding how to use this space for approximation or 
interpolation of functions. 

7.4. The Partition functions. 

Definition 7.5. We identify a function / on A with the distribution 

AeA 

where 5 V is the Dirac distribution supported at v. 

Recall that we have defined the partition function on A by 
Vx(v) := #{(ni, . . .,n N ) \ E n * a * = v -> n i G N i- 
This is then identified with the tempered distribution 
(25) T x :=J2'Px(v)6 v 
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We can then compute its Laplace transform obtaining 

(26) LT x = Y J Vx{v)e^=J{ Y ^- a 

veA aex 
We shall think of Tx as a discrete analogue of the multivariate spline Tx- 
We also have an analogue of Bx ■ Namely, setting 

Qx(v) := #{(ni, . . . ,n N ) \ 'Y^ l n i a i = v i {0, 1}}. 

(27) B x :=Y,Qx(v)S v 
with Laplace transform 

1 _ p — 2a 

(28) LB X = £ Qx(v)e- = J[(1 + e~ a ) = J[ 

veA a£X a£X 

which implies that 

(29) B x (x) = (-l) |5| 7"x(x - 2a s ) 

sex 

with a s = J2aes a - 

7.6. Some basic modules and their transforms. As in Section 12.111 
we want to consider the module Cx generated, in the space of tempered 
distribution on V, by the element Tx under the action of a suitable algebra 
of operators. 

It is convenient to choose the algebra W(A) of difference operators with 
polynomial coefficients. 

This algebra is generated by S[U], thought of as polynomials on V and 
by the translation operators r v , v G A defined by r v f(a) = f(a + v). From 
the basic formulas (@J), © we get that, under Laplace transform, a poly- 
nomial becomes a differential operator with constant coefficients while the 
translation r v becomes multiplication by e v . 

Consider the group algebra C[A] with basis the formal elements e a , a £ A. 
Define next the algebra W(V) generated by S[U], thought of as differential 
operators with constant coefficients on U and by the functions e v ,v € A. 
Additively we have W{V) = S[U] <g> C[A]. 

The algebras W(A) and W(V) are isomorphic by the isomorphism <f> de- 
fined by 

4>( T v) = e v , (f>(u) = -D u 
for v E A, u & U. Thus, given a module M over W(A) we shall denote by 
M the same space considered as a module over W(V) and think of it as a 
formal Fourier transform. 

Define S x := C[A][n aeX ( 1 ~ e"")" 1 to be the localization of C[A], ob- 
tained by inverting 5x '■= Yl a ex(^ ~ e_a ); an< ^ consider Sx as a module over 
W{V). " 
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Given a linearly independent subset A C X, let 

Aa := n « a ' n " G Z i' := (52 ra °' n ° e M+ ^' 

aeA aeA 

be the sublattice and the positive cone that they generate. Set 

£,A ■= 52 Sv 

veA A r\C(A) 

a tempered distribution with Laplace transform given by: 

( 3 °) L ^ = II TZ^r- 

aeA 

In analogy with Theorem 12.121 using [22], we get, 

Theorem 7.7. Under Laplace transform, Cx is mapped isomorphically onto 
Sx- In other words we get a canonical isomorphism of Cx with Sx as 
W(U) -modules. 

Cx is the space of tempered distributions which are linear combinations 
of polynomials times ^a, A cz X a linearly independent subset, and their 
translates under A. 

Proof. The proof is analogous to that of Theorem 12.121 We use the results 
on partial fractions explained in [22^ to show that Sx is generated by S^ 1 
as a W(V^)-module, together with the simple formula ()3U|) . □ 

Remark 7.8. There is a similar parametric case where the Laplace transform 
takes values in the algebra := C[A] [fLexC 1 ~~ e~ a ~ tla )~ 1 ]. 

Notice that the fact that the Laplace transform of Cx is an algebra means 
that Cx is closed under convolution. 

7.9. The toric arrangement. The algebra Sx has a precise geometric 
meaning. In fact C[A] is the coordinate ring of an affine variety whose 
points are the homomorphisms A — > C* called complex characters. We shall 
denote this variety by T. 

Notice that each character is obtained as follows. One takes a vector 
(f) € Uq = hom(V, C) (the complexified dual) and constructs the function 
a i — t e^l°). The elements of A* := {(f> G U~c | (4> I «) € 2iriZ} give rise to 
the trivial character. Thus T = C/c/A* is an algebraic group isomorphic to 
(C*) s . 

The expression e^' a ' has to be understood as a pairing, i.e. a map 

e (4>\a) : Tx A^C*. 

This duality expresses the fact that A is the group of algebraic characters 
of T. 

The class of a vector cj) G Uc in T will be denoted by so that the value 
of e a in is e^ a >. 
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If we fix an integral basis e% for A and set Xi := e e *, we see that C[A] = 
Cfx^ 1 , . . . jxf 1 ] is the ring of Laurent polynomials, i.e. the coordinate ring 
of algebraic functions on the standard torus (C*) s . 

The equation riaex(l — e ~") = 0> defines an hypersurface Y in T, the 
union of the kernels of all the characters e a , a G X. The algebra Sx is clearly 
the coordinate ring of the complement of Y in T, itself an affine variety. 

Given a basis b extracted from X, consider the lattice Aj, C A that it 
generates in A. 

We have that A/A^ is a finite group of order [A : A&] = | det (6) | . 
Its character group is the finite subgroup T(b) of T which is the intersec- 
tion of the kernels of the functions e a as a G b. 
We now define the points of the arrangement 

(31) P(X) := U b _ eB(x) T(b). 

For any point in P(X) we choose once and for all a representative 4> G U 
so that the given point equals e^. We will denote by P(X) the corresponding 
set of representatives. We now set: 

(32) ^:={a€l|e^ |a) = l}. 

The points of the arrangement form the zero dimensional pieces of the 
entire toric arrangement. By this we mean the finite set formed of all con- 
nected components of all the intersections of the hyperfurfaces of equations 
1 — e~ a = 0, a G X. For one such intersection the connected component 
through 1 is a subtorus T", these other components are cosets of T'. 

7.10. Sx as a module. As in Section 12.111 and following |22] for every 
k < s, let us consider the submodule Sx,k C Sx spanned by the elements 

/ 

FJ.(1 - e -a^hi 

such that the vectors Oj which appear in the denominator with positive 
exponent span a subspace of dimension < k. 

Given a point G P(X) and a basis b := {ai, . . . , a s } extracted from 
Xgj,, we have seen that belongs to the finite group T(b), in duality with 
A/A;,. The finite dimensional group algebra C[A/A&] is identified with the 
functions on T(b). 

We choose a set R b of representatives in A for the cosets A/A^ by taking 
for each coset the unique representative of the form 

^Pad, < p a < I, p a G Q. 
agfe 

Character theory tells us that the element e(4>) := | det(6)| _1 ^AeK b e ~^ A ' e ' 

has the property of taking the value 1 on the point and on all the other 
points of T(b). Moreover: 
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Proposition 7.11. e(cj))Y\ a£b (l — e a ) 1 is the Laplace transform of the 
distribution 

veC(b)nA 

Proof. By our choice of representatives, if v € C{b) n A we can write v = 
X + J2a€b n a.a, n a € N, A G Rb 



□ 



In Tlx '■= Sx / Sx,s-i the class 



WM := e6 (l-e- a ) ^' m ° d ( - Sx ' s ~ 1 ^ 

is clearly independent of the chosen representatives and is an eigenvector for 
A of eigenvalue e?. Indeed 

, ^ |det(6)|- 1 2^ ApR e -^+Ml^> e A+M 

*M = e<At " ] [ nt>(i-e- a ) ] = M 

by the independence from the chosen representatives. One can show |22j : 

(1) The W(V) module IIx is semisimple of finite length. 

(2) The isotypic components of Hx are indexed by the points of the 
arrangement. 

(3) Given € P{X), the corresponding isotypic component IIx(^) 
is the direct sum of the irreducible modules n& ^ generated by the 
classes ujf,^ and indexed by the no broken bases extracted from X^. 

(4) As a module over the ring S[U] of differential operators on V with 
constant coefficients, 11^ is free of rank 1 with generator the class 

Although we do not want to reproduce here all the details of [221; ^ is not 
too difficult to reconstruct all these statements, in particular by using the 
ideas of section 19.181 

In fact, since W(V) = S[U] <g> C[A], we have: 

n&,0 = w{v) ® C [a] ub,<t> = s[u]ujb,<i> 

Theorem 17. 71 tells us that we can transport our filtration on Sx to one on 
Cx- By formula 1|30|). in filtration degree < k we have those distributions 
which are supported in a finite number of translates of the sets A n C {A) 
where A spans a lattice of rank < k. 

Remark 7.12. In ([22]) there is a similar explicit description of the piece 
Sx,k/ Sx,k-i, for each k > 0, using the varieties of the arrangement, of 
codimension k. The proofs are in the same lines as the ones we developed 
in the simpler case of hyperplanes with Proposition 13.41 
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7.13. Local expression for Tx- Let us now consider the element vx class 
of the generating function n^ex^ ~~ e x )~ l in Tlx- Decompose it uniquely 
as a sum of elements vxj, in Tlx(4>)- By what we have seen in the previous 
section, each one of these elements is expressed as a sum 

Vx <p = E %<t> u k,<l> 

for suitable polynomials q^. Thus: 

(33) v x = ^2 E ^M^M 

<t>eP{X)beNB Xlt> I ccc(fe) 

We are now ready to state but, not yet fully justify, the main formula 
that one can effectively use for computing the partition function Px- 

Theorem 7.14. Given a point x in the closure of a big cell c we have 

(34) P x (x)= e<j> E 

<f>£P(X) beNBx^ I cCC(fe) 

Proof. Using Proposition 17.111 and formulas and (|33jl . we deduce that 
the right handside of (j34j) coincides with the partition function on the cone 
C(X) minus possibly a finite number of translates of lower dimensional 
cones. 

In particular we deduce that the partition function Px coincides locally 
with a quasi polynomial for the lattice A on the cone C(X) minus possibly 
a finite number of translates of lower dimensional cones. Also, since to get 
this result we have only used the fact that each vector in X lies in A, we can 
substitute A with the larger lattice A/n and we also get that the partition 
function also coincides locally with a quasi polynomial for the lattice A/n 
on the cone C(X) minus possibly a finite number of translates of lower 
dimensional cones. By the continuity of quasipolynomials we deduce that 
our Theorem will follow from the following Proposition which will be proved 
in section 

Proposition 7.15. On the closure of each big cell c, the partition function 
Px coincides with a quasi polynomial for some lattice A/n. 

8. Partial fractions 

8.1. Algebraic identities. We follow the approach of Szenes and Vergne, 
used in |37j (and also in [22]) to prove Proposition 17.151 and hence finish the 
proof of Theorem 17. 141 

We do not need a very fine analysis, let us consider the set X formed by 
all the positive rational multiples of the vectors in X. 

Observe that all the cones generated by subsets of elements of X coincide 
with the ones generated by subsets of X. In other words the geometry of 
C(X) does not depend on the denominators that we will introduce. 
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The following lemma is a simple computation 

Lemma 8.2. Let b be a set of linearly independent vectors in a lattice M , 
e^ a character of the lattice spanned by b, n c a positive integer for all c £ b. 
If we expand 

cGfe y eMnc(b) 
we have for y = Ylceb ^cC, k c > for all c £ b 



Notice that this lemma give immediately Proposition 17.151 in the special 
case in which X consists of the elements of a basis (contained in the lattice 
A) each repeated any number of times. 

Our next task is to reduce to this case. For this we need to show 

Proposition 8.3. Let 7 = Y2aex raCL w ^ r a £ Q; < r a < 1, then the 
function e~ 7 IIaex(l — e"^~ a )~ 1 , with e^ a character, can be written as a 
linear combination with constant coefficients of elements of the form 

JJ (i _ e (4>\c)-cyn c 

ceb 

with b a linearly independent set of elements in X. Furthermore if e^ is of 
finite order, each of the e^ has finite order on the lattice spanned by b. 

Before giving the proof of this Proposition we need to show some simple 
identities. 

Lemma 8.4. 

n 

(35) 7^ = 1^ C:=e 2 ™/". 
1 — x n 1 — fx 

i=0 s 

Proof. Take an auxiliary variable t and the logarithmic derivative (relative 
to t) 

ra +n-l n ~ 1 
T^^dt = dlog(i" - x n ) = d]og(Y[(t - C^)) 
x i=0 
n— 1 ji—l ^ 

= J2 dl °s(t - Cx)) = E 77^7iZ) dt - 
i=0 i=o 1 15 } 

Next set t = 1 in the coefficients of dt. □ 

Lemma 8.5. 

(36) i-n*= e ik- 1 ) 1 ^^ 1 -^) 

i=l 0C/C{l,...,r} i£l 
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Proof. The proof is by induction on r. The case r = 1 is clear. In general, 
using the inductive hypothesis, we have 

E n(- 1 ) |/|+1 ( 1 -^)= E EK- 1 ) 1 "-^ 1 -*)- 

0C7c{l,...,r} iel 0C7c{l,...,r-l} iel 

e ri(- i ) |7|+i ( i -^)( 1 -^) = i-n^+a-a-n^ci-^) = 

7C{l,...,r-l} iel i=l i=l 

r 

-»-n 



i=i 



□ 



A variation of this formula is the following: 
Lemma 8.6. Set t = flti z i Ui=h+i z ^ > a G C * 

!-*= e ri(- i ) |7|+i ( i -^)-«" 1 e n(- i ) |/|+i ( i -^)+ 

0C7c{l,...,M i6J 0CJc{h+l,...,r} J6-T 

(37) a- 1 E II(- 1 ) m+1 ( 1 -^)( 1 - a *) 

0C/c{h+l,...,r} iel 



Proof. This is immediate from the previous Lemma once we remark that 

i-i = i-nL«i-«(i-nUi4 D 

Prom this we get, 
Lemma 8.7. Set t = \^ =1 flLh+i ™^ < /i < r. T/ien 



i ^ ( -i)m+i 



nua-^) ^ } (i-t)n^(i-^) 

(_l)W+i (-l)l'H- 1 



(38) "^cSw/iwa-*) d-oiwd-^)'- 

If a £ C* and a ^ 1 

n - oil 

0c/c{i,...,M 

(39) 

(_l)|/|+i (-l)l'H- 1 



1 o f V (~1) |/|+1 

nr=i(i -*)(!-«*) °- i (JCJ c^r ^-^iwi-* 



-— ( y 
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Proof. The first relation follows from (|H7jl dividing by (1 — t)(l — z\)(l — 
Z2) • ■ ■ (1 — z r ) and taking a = 1. 
The second writing 

1 = _2_(1_*)- J_(l_a*) 
a — 1 a — 1 

and then dividing by (1 — at)(l — zi)(l — 22) •■■(! — z r ). □ 

Lemma 8.8. Given elements a- L E X, b = Y^i( n i/ n ) a i < n j < n > an d a- 
character e v of the lattice spanned by the a,, we can write the element 

Hf^l - e^M-^) ~ f_\ (1 - e {nWi)-a^ 

as a linear combination with constant coefficients of elements of the form 

1 

(1 _ g(C|ci)-Cl) . . . (1 _ e (C|Cr)-C r ) 

where the Ci £ X and e^ is a character of the lattice spanned by the Cj. 
Furthermore if e v is of finite order, also each e^ is of finite order. 

Proof. Start from the identity 



-Ha) _ e -Ha) H _ & {r]\a)-a\ e -(vW) 



e -(v\a) 



1 — e {r)\a)-a 1 _ e {r)\a)-a \ _ e (r)\a)-a 

from which we get that, if k < n and e^,i = 1, . . . , n are the n th roots of e 71 : 

g—ka/n k g— a/n n ^ 



_rr ^_ rr _ 

l — e {r)\a)—a 1 1 ]_ _ g (Ci |a)— a/n 11 J _ e {Ci|a)— a/n 
i=l i=fc+l 

can be expanded as a linear combination of products of elements (1 — 
e (Cil«)-o/ n ) _1 , "We now apply this procedure to each term of our product. 
The fact that if e v is of finite order also the are of finite order follows 
immediately from their definition. □ 

8.9. The main expansions. We can now give 

Proof, (of Proposition From Lemma 18.81 we deduce that we can assume 
that 7 = 0. 

Also there is nothing to prove unless we have a linear dependency nao = 
Y^i n i a ii n i € Z which we write as 

lilt"*! i \ n \ij^i n 3\ 

We substitute now ao with 6o := a o/| IL n «l an d °« with h := a.^ / 1 tt. XT j m 71 j\- 
Next apply lemma (|8.4|) . using suitable characters e**. We get a new expan- 
sion in which the b% appear in place of the aj. 

For each product Yu=o 0- ~ e^ 6 *) -6 *) -1 we have two possibilities. If the 
equations 1 — e^ bi ' b * = are incompatible we are in position to apply the 
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formula (|3ll|) to X{ = de b \i > 1 and to t = (oe b °. If they are compatible 
we apply formula (|38|) . In both cases we substitute the product with a 
linear combination of products with less terms in the denominator. We then 
proceed by induction. □ 
We can now finish the Proof of Proposition 17.151 Using Proposition 18.31 
and Lemma l8~2l we get that for any 7 = J2 a ex r « a with r a G Q, < r a < 1, 
the function e~ 7 riaex(l — e ~ a )~ 1 is a quasipolynomial in the interior of the 
big cells. From this we deduce that 

1 e~ 7 

is a quasipolynomial on the big cells translated by 7. Since as 7 varies these 
translates cover each big cell, our claim follows. □ 

8.10. E— splines. In our setting we can also consider the Euler-Maclaurin 
sums (cf. 0>[ini): 

P(v) := {(m, 

Vx,iJ,{y) must be understood as the distribution X^gAuCpf) "Px,h(v)5 v 
supported at the points of the lattice A U C(X). 
Its Laplace transform is 

v£AuC{X) a&X 

This leads to a similar theory in which now in general there are more points 
in the arrangement, in the generic case S(X) (see definition I44|) . distinct 
points, and correspondingly more exponential functions appear in the final 
formulas. We leave the details to the reader and also refer to section 19.91 
where we approach the same topic from a different perspective. 

A particularly interesting case is when the numbers e^ a = v a = x( a ) f° r a 
character \ °f A, i n particular for a character of finite order which is 1 on a 
sublattice Ao of finite index. Then we have that, if ^2 a£ x n a a = v we have 
e -T,aex n °.Va = x(v). In this case we denote Vx,n = Px,x- Thus we have 
that _ 

Px x = Yl x(v)Vx(v)5 v 
veAuc(x) 

and 

^k = iA/A r 1 y, Vx >x 

xeA/A^" 

where Vx\a is the restriction of Vx to Ao. In a similar way, by applying 
translations we can restrict to the other cosets. 



- 2^ 



- E «" 

(ni,...,njv)eP(u) 

• • • 1 n N) 1 y2 Uiai = Vi Hi<E 



10 



C. DE CONCINI AND C. PROCESI 



9. Difference equations 

9.1. Difference operators. Let us consider the space of complex valued 
functions / on the set A. We identify this space with the algebraic dual 
C[A]* of C[A] by the formula: 

(f\e a ) :=f(a). 

For v G A we define the difference operator V„, acting on functions / G C[A]* 
as: 

V„/(«) := f{u)- f(u-v). 

Parallel to the study of D(X), we start by studying the difference equations 
Vy/ = 0, where Vy := n^ey V„ as y G 8(X) runs over the cocircuits. 
Let us denote the space of solutions by: 

(41) V(X) := {/ : A - C, | Vy(/) =0,VFe f (X)}. 

As we shall see, this space is not only a formal construct, but it plays an 
essential role in the theory of partition functions. 

We want to reformulate the fact that a function is a solution of a system 
of difference equations as the property for such a function to vanish on an 
appropriate ideal Jx of C[A]. 

Notice that, the ideal I\ of functions in C[A] vanishing at 1 G T has as 
linear basis the elements 1 — e~ a , a G A, a ^ 0. If one takes another point 

the ideal 1$ of functions in C[A] vanishing at e^, has as linear basis the 
elements 1 - e - a +M4>> ; a G A, a ^ 0. 

In fact, for every a, x G A we have: 

(V a /|e*) = (/|(l-e- a )0. 

Thus the difference operator V a is the dual of the multiplication operator 
by 1 — e~ a . In this setting, we get a statement analogous to that of Theorem 
15. Ml on differential equations: 

Theorem 9.2. A function f on A satisfies the difference equation 

p(V ai , . . . , Vat)/ = with p a polynomial if and only if thought of as 
element of the dual o/C[A], / vanishes on the ideal o/C[A] generated by the 
element p(l — e~ ai , . . . , 1 — e~ ak ) . 

We have also the twisted difference operators V* G T defined by 
(Vtf)(x) := f(x) - e^^fix - a), dual to multiplication by 1 - e" a+ < a !</>>. 

We need some simple algebraic considerations on subschemes of T sup- 
ported at a point G T. Let S =: C[A]/J be the (finite dimensional) 
coordinate ring of such a scheme. 

In this slightly more general setting, we can repeat the discussion made 
in section Tl. 131 Given an element a G A, we have that e~^ a ' ^ e a — 1 is in (f>. 
Thus, in the ring S, the class of e~^ a ' ^ e a — 1 is nilpotent. The power series 
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of log(l + t) = t - t/2 + t 2 /3 - ... can be computed on t = e~( a I ^ e a - 1 
and we defines a := log(l + (e~< a I *>e a - 1)) + (a | 0). 
We clearly have: 

Lemma 9.3. The map i : a — > a is a linear map. 
a := a — (a\4>) is nilpotent. 
e- = exp(a)e^ a l^ G 5 egua/s t/ie c/ass o/e a . 

i extends to a linear map i : V — > S and then to a surjective homomor- 
phism i : S[V] — > S 1 . 

Proof. The class e -( a l0}+ a = exp(a) = ^2/, >0 a k /kl reduces to a finite sum 
lying in the image of the homomorphism i which is therefore surjective. □ 

We have thus identified S with a quotient S[V]/I by an ideal of finite 
codimension in S[V}. 

The fact that, for each a G A we have a — (a | <f>) is nilpotent, means that 
/ defines a unique point in which of course a = (a \ <j>) . That is / defines the 
point 4> (notice that the choice of 4> is not unique). 

We thus can identify, using this algebraic logarithm, the given scheme as 
a subscheme in the tangent space. This proves the following: 

Proposition 9.4. Let J C be an ideal such that C[A]/J defines (set the- 
oretically) the point e^. Under the logarithm isomorphism C[A]/J becomes 
isomorphic to a ring S[V]/I defining the point <f>. 

In particular we get a canonical isomorphism of the space of solutions 
of the difference equations given by J with the space of solutions of the 
differential equations given by I. 

The explicit formula is, given f G S* = (S[V]/I)*: 

(42) /(a):=(/| e a ) = e < a l^(/|exp(a)}. 

where (f \ exp(a)) is a polynomial in a. 

A special role is played by the points of finite order m, i.e. characters 
on A whose kernel is a sublattice A^ of index m (of course this implies 
m<j) G A*). 

As we have seen, a function f(x) appearing in the dual of C[A]/J, when 
J defines (set theoretically) the finite order point e^, is of the form f(x) := 
e ( x \'t>) g(x) with g(x) a polynomial. 

Since is of finite order, we know e^ x ' ^ takes constant values (m— th 
roots of 1) on the m cosets of the sublattice A^. Thus f(x) is a polynomial 
only on each such coset. This is a typical example of what is called a periodic 
polynomial or quasi-polynomial. 

Formally: 

Definition 9.5. A function / on a lattice A is a periodic polynomial, if 
there exists a sublattice A of finite index in A, such that /, restricted to 
each coset of A in A, is (the restriction of) a polynomial. 
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We deduce the general: 

Theorem 9.6. Let J C be an ideal such that C[A]/J is finite dimensional 
and defines (set theoretically) the points e^ 1 , . . . , e^ k , all of finite order. 

Then, the space (C[A]/J)* of solutions of the difference equations asso- 
ciated to J is a direct sum of spaces of periodic polynomials e^ a ' ^«p(a) for 
the points 4>i, each invariant under translations under A. 

Although we will not need it, for completeness we now want to prove a 
converse. Let us thus take a finite dimensional vector space Q spanned by 
periodic polynomials e^ a ' ^p(a) (for some points e^* of finite order) invariant 
under translations under A. We need the following: 

Lemma 9.7. Any non zero subspace M of Q, invariant under translations 
under A, contains one of the functions e^ a '^. 

Proof. Indeed, a non zero subspace invariant under translations under A con- 
tains a non zero common A-eigenvector. This eigenvector must necessarily 
be a multiple of one ofthee< a l^. □ 
Let J be the ideal of the difference equations satisfied by Q. 

Theorem 9.8. J is of finite codimension and Q is the space of all solutions 
of J. 

Proof. We can easily reduce to the case of a unique point (p. In fact the 
commuting set of linear operators A on Q has a canonical Fitting decompo- 
sition relative to the generalized eigenvalues which are characters of A. Our 
hypotheses imply that these points are of finite order. 

For any p € Q and b € A, we set r^p = p(x + b). Consider the function 
j(p) on A given by 

(43) j(p)(b)=p(b) = (r b p)(0). 

By definition j(T a p)(b) = j(p)(a + b) so j : Q — > C[A]* is equivariant. We 
have j'(e^ a '^)(6) = e^' 9 ^. In particular e^'l^ is not in the kernel of j. 
Since j is equivariant, Lemma 19 . 71 implies that j is injective, and we deduce 
a surjective map C[A] — > Q* . 

We give to Q* an algebra structure by defining the product of two linear 
forms as follows. For /(a) = e^ a ^^p(a) € Q we have that, by translation 

invariance, f(a+b) = E f e<°l*>j>i(a)e< 6 l*>«(&) with e< a I ^Pi(a), e <b 1 ^qi(b) 6 
Q. This can be thought of as a coalgebra structure A : Q — > Q ® Q on Q. 

By duality this induces an algebra structure for which j* is a surjective 
algebra homomorphism. Thus Q* is naturally isomorphic to C[A]/J for 
some ideal J. 

Now consider the dual map it : (C[A]/J)* — > Q. By the first theorem 
(C[A]/J)* is identified to a space of periodic polynomials. This identification 
coincides with the one given by 7r, due to the definition of j. □ 
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9.9. The difference Theorem. We go back to the equations Vy/ = 0. 
As in Section 14. 3( we can generalize our problem by considering, for any 
sequence u = {u v \ v € X} of non zero complex numbers, the products 
Ny{k) '■= rLey(l ~~ v v e~ v ) as Y runs over the cocircuits, and the ideal 
Jxiii) that they generate. We set Jx '■= «/x(L •••,!)• 

The space of functions V(X), which we want to describe, is thus identified 
to the dual of C[A]/J X . 

It turns out that also this space is finite dimensional and we are going to 
show that its dimension is the following weighted analogue 5(X) of d(X). 

We use the notations of section 17.91 Given a basis b extracted from X 
consider the lattice C A that it generates in A, A/A;, is a finite group of 
order [A : Aj,] = |det(6)|, its character group is the finite subgroup T(b) of 
T which is the intersection fo the kernels of the characters e a as a £ b. 

We now define, 

(44) 8(X) := Idet(6)|. 

We have: 
Proposition 9.10. 

(45) Yl d(X (j> )=8(X). 

c*eP(x) 

Proof. S(X) counts the number of pairs e^,b such that G T(b), or equiv- 
alently bcX^. □ 

Lemma 9.11. The variety defined by the ideal Jx is P{X). 

Proof. The proof is identical to that of lemma 14.41 and of the first part of 
Theorem l4.6| the only difference being the fact, that when we extract a basis 
b from X the equations e a — 1 = as a € i define T(b). By its definition 
<GDJ, P(X) = UT(6). □ 
Now by elementary commutative algebra we know that 

C[A]/J x = e e , eP(x) C[A]/J x (0) 
where C[A]/Jx(</>) is its localization at e^. 

Theorem 9.12. Under the logarithm isomorphism C[A]/ Jx(<ft) becomes iso- 
morphic to the ring Ax^ = S[V]/I Xij> . Thus: 

(46) C[A]/J X = ® e « mx) A Xr 
In particular dim(C[A]/Jx) = o~(X). 

Proof. Let us see what happens to the equations Ny = IlueyU — e ~ v ) as 
Y runs over the cocircuits. 

When we localize at the elements 1 — e~ v where v ^ X<j> become 
invertible and hence can be dropped from the equations. As for 1 — e~ v 
where v € X^ we have that 1 - e~ v = v(l - ^ fc >i (-v) k /( k + 
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Obviously (1 — Ylk>i (~v) k /{k + 1)!) is invertible hence we can replace 
the equation Ny by My = Y\ ve y v. We obtain the equations defining A X(j> 
and this proves the claim completely. 

□ 

It is convenient to single out, in this direct sum, the term relative to 
e<t> = 1 which is D(X) and write 

(47) E(X):=(B e , £P{xle ^D(X <)> ), V(X) = D(X) © E{X). 

9.13. The parametric case. One could treat in a similar way the para- 
metric case. 

We take a sequence y_ := {v a \ a € X}. 

The main difference is that now the points of the arrangement are defined 
as follows. 

Given a basis b extracted from X instead of the finite subgroup T(b), 
intersection fo the kernels of the characters e a as a £ b we have to consider 
the set T„(b) where v a e~ a — 1 = this is a coset of the subgroup T(6), still 
consisting of | det(6)| elements. 

For a generic sequence u_ the ideal Jx{m) defines S(X) distinct and reduced 
points, while for special values we may have less points but in each such point 
a similar space of quasi periodic polynomials. 

So now we define (as for (|3*T|0 

Pu(X) := U b _ eB(x) T„(b). 
For any G Pv{X) set: 

(48) I^):={a6l|e H *=4 

We have as usual the analogues of the results in the non parametric case: 
Lemma 9.14. The variety defined by the ideal Jx(e) is Pu(X). 

C[A]/J x (u) = ® e ^ PAX) C[A]/J x (v)^) 

where [C[A]/J x (e)](4') 1S its localization at e^. 
Finally 

Theorem 9.15. Under the logarithm isomorphism C[A]/J x ((j>) becomes iso- 
morphic to the ring S[V]/I X<I> . Thus: 

(49) C[A]/J X ( E ) - ® e * ePv _ {x) S[V)/I x ^v). 
For all v_, we have dim(C [A] /Jx(k)) = 

Of course also in this case we could perform the polygraph construction, 
now the base would be a torus C . We still have a flat family and the 
polygraph is a Cohen Macaulay variety as in Remark (|4.7|l . 
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One can use the parametric case as follows. Suppose that we pass, from 
the lattice A to a finer lattice A/n for some positive integer n then, for a G A 
we write in C[A/n] 

l-ue~ a = J] (l- 7 e- a /«) 

Thus if we take our list X = {a\, . . . , azv} with parameters u = {u\, . . . , z/jy} 
in A and consider it in A/n, we get an ideal J^{u) C C[A/n]. The same 
ideal is also associated to the list 

X^ n ' = {ai/n, . . . , ai/n, ax /re, . . . ,aj^/n} with parameters 

V v ' V „ ' 

n n 

v {n) := {7i,e7i,...,e n ~ 1 7i,---,77V,---,e n ~ 1 7iv} 

where, for each i = 1, . . . N, ji is an n-th root of V{ and e = exp(27ri/re). 

One has to remark that A/re is the character group of a torus T l / n which 
maps surjectively it : T l / n — > T to T, with kernel -fC n the character group of 
(A/re)/A. 

The points of the arrangement associated to X in T 1//n are tt 1 (P(X)) 
and a union of cosets of K n . 

The algebra C[A/n]/J x (v)C[A/n] = C[A/n]/ J x{n) {v {n) ). Thus even if we 
start from the trivial parameters = 1, once we pass to a finer lattice we find 
as parameters the roots of 1. 

9.16. A realization of C[A]/Jx- As for the box spline the functions in 
D(X) play a basic role, so do the functions in V(X) in the theory of linear 
diophantine equations. In fact they describe combinatorially the partition 
function Vxiy) associated to X. 

We go back to the space Tlx '■= S x / S x>n -i of polar parts. This we have 
decomposed, as Tlx = ffi<^eP(x)-Px'j > (0)i through the points cj) £ P(X) into 
local modules Px^ifi) of polar parts for the affine arrangement centered at 
0- 

The element vx class of the generating function HaexC^ ~ e a ) m ^ x 
decomposes into a sum of local elements vx^ ■ 

We consider next the submodule Qx generated in Tlx by vx under C[A] 
and we deduce, reducing to Theorem (|6.3I) and using Theorem 19 . 1 21 that : 

Theorem 9.17. The annihilator of vx is the ideal Jx generated by the 
elements Ny = Yla^yO- ~ e_a )> as Y runs over the cocircuits. Thus Qx = 
C[A}/J X . 
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We obtain a canonical commutative diagram, made of isomorphisms, 
compatible with all the identifications made: 

C[A]/J X — =-> ®^P(x)A x , 



9.18. Prom volumes to partition functions. We apply the previous the- 
ory to the partition function Vx( a )- 

Consider the class of a product riaex(l ~~ e"")^ 1 in Tlx- For each point 
4> of the arrangement we project to the corresponding isotypic component 
Ux(4>) which is generated by the elements w&,</>- 

Since u>b,<p is an eigenvector for A of eigenvalue e^, using the commutation 
relations with the derivatives, we clearly have that, for all a <G A, the element 
e° — e^' a ^ acts locally nilpotently on the modules Uxirf') • It follows that 
the action of e a on IIx(</>) is given by an operator of the form e^ a ' ^u a with 
u a = e a with a locally nilpotent. 

The map a i— > a + (a \ 4>) is linear. It induces a linear map of V = A C 
into linear operators on Jlx(<j>) whose image consists of mutually commuting 
operators. In this way we obtain a S[V] module structure on Ux((f>)- The 
commutation relations with the operators of S[U] then give us an action of 
the usual Weyl algebra W(U). 

Proposition 9.19. As a W(U) module Hx(4>) * s isomorphic to the module 
of polar parts for the list X<p with parameters (ft := {<p a = —(<p\a)} for each 
a € X<f, by an isomorphism characterized by 

3<t>(vb,4>) = Ub 

for each no broken basis extracted from . 

Proof. The fact that the class ujb,4> is an eigenvector is an eigenvector for A 
of eigenvalue and the definition of the S[V] module structure imply that 
LUb,<p is an eigenverctor for V of eigenvalue <p. Thus we get a surjective map 
of W(U) modules 

7*:i%(£)-(n*)* 

with 70 = ^b,4> for each no broken basis in X^. Since both modules are 
free 5[C/]-modules and 7^ takes a basis for Px^{4>) to a basis for (Hx)<j> we 
deduce that it is an isomorphism, is the inverse of 7^. □ 

By the previous discussion e a acts on ILx(<j>) by the operator ' a ^e a with 
a locally nilpotent. We distinguish two cases. 

If & {<t>\a) _t lj then (i_ e -ayl 

gives an invertible operator in II x (</>)• 

If on the other hand, e^' a ^ = 1, the operator a/(l — e~ a ) is invertible 
(its power series is written in term of Bernoulli numbers). Thus the image 
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under of the component of riaexU — e a ) 1 equals 

n i-e-Uia n T^sb-vwr^Q* n^-^i")) -1 

where 

^ = II 1 _ e -a-{4>\a) II 1 _ e -a 
a^X^ aeXj, 

is an invertible operator, on Px (<ft), expressed as a power series in the oper- 
ators a and hence locally given by a polynomial. We deduce 

Proposition 9.20. Under the isomorphism j of coordinates j^: 

j ■ ®<t>£P(X)H-x(<t>) -> ©06P(X)^(0) 

we have the transformation 

(50) j: J] (l _ e-)- 1 -» Y, 0* II^-^I ))" 1 - 

a€X 4>£P(X) aex^ 

We now apply Laplace transform and deduce the final general formula 
expressing the partition function as a sum of transforms of the local mul- 
tivariate splines. In order to justify our results let us make some remarks. 
We have made an identification of Ux(4>) with Px (0) as W(U) modules. 

In particular we have that Ux(4>) is a free module, over the algebra of 
differential operators with polynomial coefficients, in the classes of the el- 
ements e(0)J7 ae6 (l — e _a ) _1 while IIx(<^) is a free module, over the same 

algebra, in the classes of Ilae^ ~~ I a ))~ 1 - 

Take then M±,M2 to be the two free submodules, over the algebra of 
differential operators with polynomial coefficients, of Sx re. Rx generated 
by these elements. Define by j : M2 — > M\ to be the module isomorphism 

mapping Uaeb( a ~ (0 I to e^ILe^ 1 ~ e _a ) _1 - 

From proposition 17. Ill we see that M\ is the span of the Laplace trans- 
forms of distributions of type 

Y ptCaOldetQOr 1 y 

<l>eP{x),bex^ v<=c(b)nA 

while M2, by © and l|l(JI) . is the span of the Laplace transforms of distri- 
butions of type 

If we restrict both types of distributions on the set of regular points of 
C(X) we see that, M2 can be identified with a space of smooth functions, 
which are locally linear combinations of polynomials times exponentials, 
while M2 can be identified to functions of the regular points in A which are 
locally quasi polynomials. 
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Theorem 9.21. Under these identifications the map j consists simply into 
restricting the functions to the points in A. 

Proof. Since at the level of functions the map is a linear isomoprhism com- 
patible with the multiplication by the polynomials S[V], it is enough to 
verify the statement on the generator XXeC7(6)nA 4>( v )^v and e^'^xc(6) f° r 
which it is clear. □ 

We deduce: 

Theorem 9.22. On the intersection of A with the open set of regular points 
we have: 

(51) V X = Q^x,^ 

0eP(x) 

Proof. The explicit formula is a consequence of formula (|5U|) plus the previ- 
ous discussion which implies that the two sides coincide, since both functions 
are restrictions of quasi polynomials. □ 

From this formula we can deduce one valid everywhere using the method 
of Jeffrey-Kirwan residues. 

Theorem 9.23. Given a point x in the closure of a big cell c we have 

(52) P x (x)= E Pb,X,(-x)). 

<t>eP{X) beNB X(t> I ccc(b) 

Proof. Both terms of the equality are continuous on the closure of c, they 
coincide in the interior by the previous proposition, hence they are equal. □ 

10. Some Applications 

In this short section we would like to give a streamlined presentation 
of some of the applications. All the results are taken from the papers of 
Dahmen and Micchelli, or from |18j . with minor variations of the proofs. 

For further details and more information the reader should look at the 
original literature. 

The main steps are the following. 

10.1. Discrete convolution. Let us first analyse the discrete convolution: 

Bx *d V = E Bx ( x + A )^(- A ) 
AeA 

and prove that: 

Theorem 10.2. When p G D(X) we have that also Bx *dP € D(X). 

This defines a linear isomorphism F of D(X) to itself, given explicitly by 
the invertible differential operator Fx '■= Yiaex ^~T> — ~- 
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Proof. A way to understand this convolution is by applying the Poisson sum- 
mation formula to the function of y, Bx(x + y)p(—y). Its Laplace transform 
is obtained from the Laplace transform e x Y\ ae x(^ ~ e~ a )/a of Bx(x + y) 
by applying the polynomial p(x) as differential operator. 
In our definition of Laplace transform we have 

Lf(0 = (27r)"/ 2 /(iO 

where / denotes the usual Fourier transform. So Poisson summation formula 
gives: 

£l#m) = X>(A) 

H A 

where fx runs in the dual lattice A*, of elements for which (fx | A) G 27riZ, 
VA G A. 

Thus if we are in the situation that Lcf)([i) = 0, V/x / 0, |iG A* we have 

L#0)=J>(A) 

A 

The main observation of Dahmen and Micchelli is that: 

Lemma 10.3. The Laplace transform of Bx{x + y)p(—y) vanishes at all 
points fi ^ 0, fjL € A*. 

Proof. We may assume that p{x) is homogeneous of some degree k. The 
evaluation of p(x) against flaex 1_ a can be understood as follows. Each 
factor can be expanded in power series. We select from at most k 

factors the homogeneous parts of some degrees hi > so that J2i hi = k 
and evaluate p against the resulting monomial, then we multiply by the 
remaining factors and sum over such choices. Now, if the factors we have 
chosen correspond to a cocircuit the evaluation of p on the monomial is 0. 
If instead this is not the case we have still a product J^aGA l ~a w here the 
elements A span. Thus, if fx ^ there is at least one a € A which does not 
vanish at fi. But now clearly 1 — e~ a and hence (1 — e~ a )/a vanishes at fi 
and the Lemma follows. □ 
We go back to the proof of our Theorem. We have shown that in our 
case, Poisson summation degenerates to the computation at 0. Taking a 
polynomial in derivatives and computing against a function and then at 
is just duality thus: 

p( x )( e * Yl 1^0(0) = <p|e* J] — —) = (Fxp\e x ) = F x p(x). 
a a 



Since D(X) is stable under derivatives and Fx is clearly invertible both 
our claims follow. □ 
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10.4. Paving the box. Take the box B(X) := {J2aex t a a I < t a < 1}, 
which is the support of the box spline. 

We start by giving a nice decomposition of B(X) into suitable paral- 
lelepipeds. In order to present it we need the following: 

Lemma 10.5. If a point v := Yl a eX v a a i < v a < 1 is in the boundary of 
B(X) the set A := {x | < v a < 1} does not span V. 

Proof. If A spans let us extract a basis b\, . . . , b s from A, take < e small. 
Then the set of points v + Ylt=i 1^1 < e i s an open ball contained in 
5(4 

□ 

Given A £ A and a set of linearly independent vectors b := {b±, . . . ,bh} 
from X define: 

h 

n A (6) := {X + ^Uh}, < ii < 1. 
i=i 

Proposition 10.6. B(X) can be paved with parallelepipeds, of the form 
H\(b) where b runs on the set of all bases extracted from X and X G A 
depends on b. 

Proof. Suppose that X = {Z,y} and we have paved B(Z) (by induction) so 
that its boundary is paved by some faces of these parallelepipeds. Consider 

B(Z) y := {p G B(Z) \p + ty£ B(Z), Vt > 0}. 

We have a map ir : B(Z) y x [0, 1] — > B(X), (p, t) i— > p + ty. We easily see 
that: 

(1) 7r is a homeomorphism to its image. 

(2) B(Z) n ir(B(Z) y x [0, 1]) = B(Z) y 

(3) B(X) = B(Z) U ir(B(Z) y x [0,1]) 

We are now going to pave B(Z) y by s — 1 dimensional parallelepipeds 
IT\(c) with c = {ci, . . . , c s _i} running on all sets of s — 1 linearly independent 
vectors in Z which are together with y form a basis. 

Consider the set H y of all hyperplanes H generated by subsets of Z and 
not containing y. Given such H take the unique linear form (f)jj vanishing 
on H and such that {4>h\v) = 1- Set 

X H, y x - 

x£Z\4> h (x)>0 

We claim that 

B(X) y = U H en y B(Z n H) + \ H , y 
and this is a paving. Remark that (pH takes its maximum value on B{Z) 
in the point \H,y Since if b > and u G B(Z f] H) + A^y, + &y) = 

4>H(^H,y) + b > (j>H(^H,y) w e get the inclusion B(Z n i?) + A# i2/ C B(X) y . 

To see the converse observe that by Lemma |10.51 B{X) y is a union of 
polytopes of the form B(Z n H) + \x with H G H y and fx = ^2 aeA s a a with 
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A = Z - (Z n H) and e a G {0, 1}. Fix a e A. Write y = h + ^ H {a)~ l a. 
Take v in the relative interior of B(Z n i/). We have 

■u + = (v + i/i) + (n - e a a) + (e a + t<f> H {ay 1 )a. 

If i is sufficiently small then v + th G -B(Z n -ff). Furthermore if e a = and 
<t>H{p) > or e a = 1 and (f>H(a) < 0, < e a + tcj)H(a)^ 1 < 1. Thus this 
point lies in B(Z) giving a contradiction. 

By induction, each B(ZPiH) is paved by s — 1 dimensional parallelepipeds 
II^(c) with c = {ci, . . . , c s _i} running on all bases of H extracted from ZilH. 
Since {c, y} := {ci, . . . , c s _i,y} is a basis of X and all bases containing y 
are so obtained, we get the desired paving 

n(B(Z) y x [0,1]) = U A ,JI A ({c,y}). 

From this our claim is immediate. □ 
As a simple application we obtain 

Proposition 10.7. Let xq be a regular point. Then (B(X) — xo)CiA consists 
of 5(X) points. 

It is easy to see that each parallelepiped 11^(6) translated by a regular 
point, intersects A in |det(6)| points. Summing over all parallelepipeds we 
obtain our claim. □ 



10.8. Linear independence. We now assume that we are in the unimod- 
ular case. In this case 5(X) = d(X). Choose a regular point xq consider the 
d(X) points P(x ) := (B(X) - x ) D A = { Pl , . . . ,p d{x) }. 

Proposition 10.9. Evaluation of polynomials in the points in P(xq) es- 
tablishes a linear isomorphism between D(X) and the C d ( x ^ (or if we 
restrict to real polynomials). 

Proof. Since dim(D(X)) = d(X) it suffices to prove that, a polynomial 
p(x) G D{X) vanishing on these points is identically zero. From the formula 
p(x) = J2 a e\ B(x+a)p(—a) we see that, for given x the only terms appearing 
in the sum are the ones where x + a G B(X) or a G B(X) — x. Thus for x in 
a small open neighborhood of the given point xq the only terms appearing 
are the ones in which a G {pi, . . . ,Pd(x)}- If by way of contradiction, a 
polynomial p G D(X) vanishes on these points we have that p(x) = 
proving the Proposition. □ 
Now we can prove the Theorem on linear independence of the translates 
of the box spline: 

Theorem 10.10. For X unimodular and any, non identically 0, function 
f{\) on A we have: 

Y J B(x + a)f(-a)^Q. 

a€A 
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Proof. Assume /(ao) ^ for some ao £ A. We can find then, a regular point 
xq, such that ao € "P(xo) is one of the points p% previously defined. Thus 
there is a non zero polynomial p(x) € D(X) coinciding, on V(xq) with /. 

Then J2aeA B ( x + a )f(~ a ) = p( x ) + on the set B ( x ) ~ x o- □ 
Unimodularity is a necessary condition for this theorem in fact one has: 

Proposition 10.11. Discrete convolution maps V(X) into D(X) with ker- 
nel E{X) (see Formula JJTtyj. 

Proof. The first statement follows from the following identity, for Y C X we 
have: 

D Y {B X *d a) = B x /y *d Vya 

The second by the fact that the kernel of the discrete convolution is in- 
variant under translation thus it suffices to verify it for the functions as 
(f> varies over the points of the arrangement. Let Y C X be a basis of a 
sub-lattice Ao on which <j> is 1. We have Bx *d <P = Bx/y * By *d <fi an d 
By *d 4> = since X^aeA/A ^( a ) = ^ (* ne sum * s over cose t representatives). 

□ 

We apply the previous theory to the partition function Vx( a )- 
Proposition 10.12. 

T x (x) = J2'Px{a)B x (x - a) 

aeA 

Proof. Compute the Laplace transform L(J2 ae j^ Vx{p)Bx{x — a)): 

= ^Vx(a)e-«L(B x ) = ]J ^ R ^ = R \ = LT X . 

aeA aex aeX aeX 

□ 

In the unimodular case, where we have the linear independence of trans- 
lates of Bx , we recover the results of section 19.181 
On D(X), Fx is invertible and its inverse is: 

xex e 

We have by definition QFp = p = FQp on D(X). Take a big cell c over 
which Tx coincides with some polynomial p c £ D(X). Set q c := Qp, we 
have on c 

Tx(x) = FQp, = QPc(a)B(x - a) 

aeA 

since Tx(x) = YlaeA^xi^Bxix — a) we have by linear independence: 



(53) 



QPc{a)=V x (a). 
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Notice that Q is like a Todd operator, its factors can be expanded using 
the Bernoulli numbers B n by the defining formula: 

- V — (-D ) n 

i _ e -D x ~ Z> n ^ u *> 

k=0 

Formula allows us to pass from the formula for the volume to one for the 
partition function (a special case of formula l9.23j) . 



Part 3. Residues 



11. Residues 

11.1. Cohomology. For the computations of volumes and partition func- 
tions we want to apply a cohomological method, like the usual method in 
one complex variable for computing definite integrals. 

With the notations of section I2.11| we start from the affine algebraic 
variety Ax, with coordinate ring Rx- Using the De Rham's Theorem due 
to Grothendieck, in order to compute the cohomology of Ax with complex 
coefficients we can use the algebraic de Rham complex (tt x ,d) (cf. |27j). 
Here d is the usual de Rham differential while for each < k < s, the 
algebraic differential forms of degree k are: 

Q X ■= { fii,...,i k ( x )d x ii A ■■■ Adxi k \ fi 1) ... ) i k {x) G Rx} 

l<ii<i2<-<«fc<s 

where {xi, . . . ,x s } are coordinates in U. 

Let us look more closely to what happens in degree s, the top degree. 
Define a homomorphism, uniquely determined by the choice of the volume 
form dx\ A • • • A dx s , 

by ixifdxi A • • • A dx s ) = f. We have that ixid^x' 1 )) is the space d(Rx) 
spanned by all partial derivative of elements in Rx- 

Definition 11.2. Let Hx be the space spanned by the elements Mb := 
riae6 a_1 ' as - var i es among the bases extracted from the list X. 

From the computation of the cohomology of Ax for which we refer to 
|32j . we get the following result which also follows immediately from the 
discussion of the next section and the expansion in partial fractions. 

Theorem 11.3. (1) We have the decomposition in direct sum: 

Rx = H x ® d(Rx) 

In particular, using ix, the space Hx can be identified to the space 

ir'iAx). 

(2) The elements Mj,, as b varies in the set MB{X) of no broken bases 
are a basis of Hx- 
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If we take any basis b = {b±, . . . b s } of V we have 

dbi A • • • A db s = det(b)dxi A • • • A dx s 
(this is gives the definition of det(b)). It is then clear that 
s 1 

det(fe) r-dxx A • • • A dx s = d\og(b\) A • • • A d\og(b s ). 

■ 1 ®i 

Using ix we deduce: 

Proposition 11.4. The cohomology classes of the forms 

u>b '■= dlog(ai(x)) A • • • A dlog(a n (x)) = dlog(a\) A • • • A dlog(a n ) 

as a\, . . . ,a n varies in NB{X), form a basis of the cohomology H n {Ax)- 

11.5. Local residue. From the expansion in partial fractions follows eas- 
ily that Rx = Hx + d(Rx), thus in fact Theorem 111.31 is a consequence 
of Proposition 111.41 This can be proved directly using the following local 
computation. 

Let b = {bi, . . . , b s } be a basis extracted from X. To b we associate an 
injection 

jb : Rx -> C[[ui, . . . ,u s }][(ui ■ • - lis) -1 ] 

(C[[«i, . . . , u s ]] [(ui ■ ■ ■ Ug)^ 1 } is the ring of formal Laurent series in the vari- 
ables ui, . . . , u s ), defined by 

jb(f(h,...,b s )) = /(ui,uiu 2 ,...,uiu 2 ---n s ). 

This is well defined since given a € X, if k = r jb( a ) is the maximum index 
such that a G . . . , b s ), we have 

s s j 

a = ^2 a i b j = u i " " " u k(ctk + ^2 OLj Yi u i) 

j=k j=k+l i=k+l 

with ah 7^ so that a" 1 can be expanded as a Laurent series. Clearly this 
map extends at the level of differential forms. 

Given a top differential form u G we define resb(io) as the coefficient 
of (u\ ■ ■ ■ u s )~ 1 du\ A • • • A u s in jb(u)- 

The main technical result from which all the important results follow is 
given by 

Lemma 11.6. Let b and c be two no broken bases extracted from X then 

resb(uJc) = h,c- 

Proof. First notice that jb(^b) = dlog ui A ■ ■ • Ad log u s so that resbiyJb) = 1. 

Now notice that if c ^ b there are two distinct elements c,c' G c such that 
lb(c) = 7b(c'). Indeed since in a non broken basis the fist element is always 
also the first element in the list X, 7b(ci) = 1. If {c%, ■ ■ ■ c s ) ^ (62, • • • b s ) 
then there exists an index i > 1 with 75 (cj) = 1. If (02, • • • c s ) = (62, . . . b s ) 
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then both c — {ci} and b — {bi} are no broken bases of (02, . . . c s ) extracted 
from Ifl (c2, . . . c s ) and everything follows by induction. 

We than write each c in the form ]X =1 Uif with /(0) 7^ so that (ilog c = 
^l°g(Ili=i u i) +dlogf. Expanding the product j&(wc) we get a linear com- 
bination of forms, all terms containing only factors of type dlog(Y\i = i u i) 
vanish since two elements are repeated, the others are a product of a closed 
form by a form dlog(/) with /(0) 7^ which is exact, so all these terms are 
exact and the residue is 0. 

□ 

As a first consequence we deduce: 

Theorem 11.7. The cohomology classes of the forms LOb, as b varies among 
the no broken bases are linearly independent. 

For any top differential form ip, denote by [ip] its cohomology class we 
have: 

(54) [ip] = resb(ip)LUb 

b£AfB{X) 

Proof. Since clearly the map jb takes exact forms to exact forms we deduce 
that resb factors through H s (Ax)- In view of this everything follows from 
Lemma lll,6l 

Remark 11.8. It is not difficult to verify that we get the same cohomology if, 
instead of taking forms with coefficients in Rx we take coefficients of type 
f /d k with / any function holomorphic around 0. 

We are now going to define the total residue Tres and the residues resb- 
With the use of ix these operators can be defined either on algebraic 
differential forms of degree s or on functions. 



Definition 11.9. Given / € Rx, Tres(f) is the cohomology class of the 
form fdx\ A ■ ■ ■ A dx n . 

We can now reformulate formula (|54|) as: 

Tres(f) = ^ resb(f)[uJb\. 

b£NB(X) 

11.10. Residues and Laplace transform. Recall that in formula (|15|) we 
have introduced the polynomials Pb,x, m order to give an expression of the 
multivariate spline Tx ■ 

Theorem 11.11. For every b € AfB(X), 

e (y\ x ) 

(55) Pb,x(-y) = det^resbi-j^). 

- d(x) 

with d{x) = Uaex a - 
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Proof. We begin remarking that formula (|55|) makes sense since, if we expand 

e <2/F> /d(x) with respect to the variables y = {y%, . . . ,y s }, we get a power 

series whose coefficients lie in Rx- 

In order to prove this Theorem we need some properties of Tres. The 

first property of Tres, which follows from the definition is that, given a 

df 

function / and and an index 1 < i < s, we have Tres(— — ) = 0, hence for 



' dx 



two functions /, g: 

Tres{ H 9) = ~ Tres{f ^- 
In other words for a polynomial P: 

(56) Tres{P{d x )U)g) = Tres(fP(-d x )(g)). 

We shall use the preceding relation (|56j) for the function / = e^' 3 ^ for 
which we have: 

(57) P{d x )e^ = P(y)e^ x \ 

The second simple property is that, given a basis b extracted from X and 
a function / regular at 0, we have: 

(58) Tresi s f ) = f(0)Tres( ' 



Uaeb a ( x ) ' Ua& a ( x ) ' 

We get 

p {y\x) i 

= E ^ ^U^aL) ) )= E Pk,x(-y)Tres( 1 ) 
(59) = E ddW^-^N- 

From this the theorem follows. □ 

11.12. Partition functions. The same method can be applied to partition 
functions. In ([22]) we have computed the full cohomology of the toric ar- 
rangement, nevertheless this is not strictly necessary for the residue compu- 
tations which are essentially local. In the formula (|34|) one has contributions 
localized at a point of the arrangement and for a no broken basis b C X^. 
They can be computed again using the operator resb,<j, in this case. 

This means that, we restrict a function or a form to a neighbourhood of <j>, 
use logarithmic coordinates so that the divisors 1 — e~ a = 1 which appear in 
this neighbourhood coincide with the linear hyperplanes of the arrangement 
associated to X^ and finally compute the local residue at b for this local 
hyperplane arrangement. 
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By formula (JSHJ and Proposition 19 . 2( )1 we get that the class of the Laplace 
transform riaexU ~~ e~ a )~ l in the module Hx of polar parts decomposes 
into a sum of local factors, corresponding to the summands ^x,4> an d then 
ILy,</> = ®S[U] [riaGb( a ~~ (^l a )) -1 ] under the logarithm map. With this 
isomorphism j we have 

j(II( 1 - e " a )- 1 )= £ £ Wdlfc-W ))" 1 ] 

aex ^eP(x) beMB X(t> a&b 

which bv 19,211 gives a formula for the partition function. We claim that we 
have 

e <wl*> 

(60) %x+(-v) = &t(6)req^( _ e _ a(g) _ ( ^ |a)) )- 

In order to prove this, we start by making a change of coordinates x = 
2+(y|0) so that we center the point at and we pick out the factor = e^'^ 
and then showing, using the language of proposition 19 . 21)1 that: 

Lemma 11.13. res^f) = resb,<j>{j(f)) for any f G M 2 . 

Proof. Due to the properties of Tres it is enough to prove it on the genera- 
tors where we see that, around the point (j> we have 

n f-e-«) = 1 det ® I" 1 £ me ~ x IB 1 - ^ 

= | det®!" 1 £ 0(A)e' A [J a^ 1 /i(x), /i(c/>) = 1. 

aeb 

Since by definition 0(A) = e^'^, at the point we have that 

I det®!" 1 ^ 0(A)e-< A I ^> = 1. 

This proves the claim by formula (|58|) . □ 

Now we go back to formula We have that 
(61) 

res »Ai(n ^)J > =res b Ae {y]x) %xMH(a- (01a))- 1 ]) 

Change coordinates, so to center at 0, getting x = z + and so 

re^(e^l^q^(y)[n(« " & I a »^) = ^(ST^W-v)^ 

agfe 



Summarizing These formulas, together with the local computation of 
residues form an effective algorithm to compute the functions we have been 
studying. 

It remains to discuss a last algorithmic point. 
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In order to compute the Jeffry-Kirwan residue, at a given a point p G 
C(A), it is necessary to determine a big cell c for which p£c. 

In general, the determination of the big cells is a very complex problem. 
For our computations it suffices much less. 

Let us take thus simply a point q internal to C(A) and not laying on any 
hyperplane generated by n — 1 vectors of X. This is not difficult to do, and 
let us consider the segment qp. 

This segment intersects these hyperplanes in a finite number of points, 
thus we can determine an e sufficiently small for which all the points tp + 
(1 — t)q, < t < e are regular. 

If we take one of these points go it lays in a cell for which p is in the 
closure. 

At this point, for every no broken basis, we must verify in simple way if 
qo lays or not in the cone generated by the basis. 

12. Minimal models 

12.1. Geometry of residues. For the moment our definition of residues 
resbifj is purely algebraic. In fact its true geometric meaning is based upon 
a general definition of the notion of residue in several dimensions. 

This section is quite independent of the rest of the paper and can be used 
as an introduction to the theory developed in JH] and 20 1 . 

The first point to be understood is that, the non linear coordinates Ui 
used in section lll.5| represent local coordinates around a point at infinity 
of a suitable geometric model of a completion of the variety Ax- In fact 
we are thinking of models proper over the space U D Ax in which the 
complement of Ax is a divisor with normal crossings. In this respect the 
local computation done in section fl 1.51 corresponds to a model in which all 
the subspaces of the arrangement have been blown up, but there is a subtler 
model which gives rise to a more intricate combinatorics but possibly to 
more efficient computational algorithms, due to its minimality. 

13. Irreducibles and nested sets 

13.1. Irreducibles and decompositions. The notions that we are about 
to give are of combinatorial nature (cf. 25 ) but we develop them in the 
following context: 

As usual let us consider a list X := {ai, . . . , ax} of non zero vectors in V 
which in this section we assume to be a complex vector space. 

Given a sublist A C X the list A := X n (A) will be called the completion 
of A. Thus A is complete if A = A. 

The space of vectors (f) G U such that {a\4>) = for every a € A will be 
denoted by A-^. Notice that clearly A equals to the list of vectors a € X 
which vanish on A . 
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From this we see that we get a bijection between the complete sublists of 
X and subspaces of the arrangement defined by X. 
We give the main 

Definition 13.2. Given a complete set A C X, a decomposition is a 

decomposition A = A\ U A 2 in non empty sets, such that: 

{A) = {At) © (A 2 ). 

Clearly the two sets A±,A 2 are necessarily complete. 

We shall say that a complete set A is irreducible if it does not have a non 
trivial decomposition. 

If A = A\ U A 2 is a decomposition of a complete set and B C A is 
complete we have B = B\U B 2 , where B\ = A\ PI B, B 2 = A 2 n B. Also 
(B) = (Bx) (_B 2 ), and we have 

Lemma 13.3. B = B±U B 2 is a decomposition, unless one of the two sets 
is empty. 

We deduce immediately: 

Proposition 13.4. If A = A\ U A 2 is a decomposition and B C A is 
irreducible, then B C A± or B C A 2 

From this get: 

Theorem 13.5. Every set A can be decomposed as A = A± U A 2 U • • • U Ak 

with the Ai irreducible and: 

(A) = {A 1 )®{A 2 )®---@{A k ). 

This decomposition is unique up to order. 

Proof. The existence of an irreducible decomposition follows by a simple 
induction. 

Let A = B\ U B 2 U ■ ■ ■ U B^ be a second decomposition. Proposition 113.41 
implies that every B{ is contained in an Aj and viceversa. 

Thus the Aj's and the B^s are the same up to order. □ 
A = A\ U A 2 U • • • U Ak is called the decomposition into irreducibles of A. 

Example 13.6. An interesting example is that of the configuration space 
of s-ples of point in a line (or the root system A s ^\). In this case X = 
{zi - zj\l <i<j<s}. 

In this case, irreducible sets are in bijection with subsets of {1,... ,s} 
with least 2 elements. Indeed give one such subset S it corresponds to the 
irreducible Is = {zj — zj\{i,j} C S}. 

Given a complete set C, the irreducible decomposition of C corresponds 
a sequence of disjoint subsets £1, . . . , Sk of {1, . . . , s} with least 2 elements. 
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13.7. Nested sets. We define now the basic notion of nested set. 

We say that two sets A, B are comparable if one is contained in the other. 

Definition 13.8. We shall say that a family <S of irreducibles A^ is nested 
if given elements Ai 1 , ■ ■ ■ , Ai h £ S mutually incomparable we have that 
C := A 1 U A 2 U ■■■ U Ai is complete and C := A 1 U A 2 U • • • U A t is its 
decomposition into irreducibles. 

Remark 13.9. If Ai, . . . ,A^ is nested we have that UiAi is complete. In fact 
this union can be obtained taking the maximal elements, that are necessarily 
non comparable, and then applying the definition of nested. 

We are in particular interested in maximal nested sets, which we denote 
by MNS. 

Nested sets can be inductively constructed combining the following in- 
ductive procedures. 

We can use two simple inductive ways to construct nested sets. The proof 
is left to the reader: 

(1) Suppose we are given a nested set S, an minimal element A E S and 
a nested set V whose elements are contained in A. Then we have 
that S U V is nested. 

(2) Suppose we are given a nested set S, a complete set A containing 
each element of S. Then if A = A\ U • • • U A^ is the decomposition 
of A into irreducibles, S U {j4i, . . . , A&} is nested. 

In the case of Example 113.61 one can easily verify that nested sets cor- 
respond to families of subsets in {l,...,s} each containing at least two 
elements and such that any two subsets of the family are either comparable 
or disjoint. 

Theorem 13.10. Assume that (X) = V. Let S := {A 1 , . . . ,A k } be a MNS 
in X. 

Given A € S let B\, . . . , B r be the elements of S contained properly in A, 
and maximal with this property. 

(1) C := B\ U • • • U B r is complete and decomposed by the B{. 

(2) dim(A) = dim(C) + 1. 
(8) fc = dim(F). 

Proof. (1) is the definition of nested set since the Bi, being maximal, are 
necessarily non comparable. 

(2) Let us consider (C) = ffi- = l( B i) c (^)- 
(C) ^ (A) otherwise, since C is complete, by the definition of nested set, 
we must have A = C. This is absurd since A is irreducible and the Bi's are 
properly contained in A. 

Therefore there exists an element a € A \ a ^ (C). Let us denote by 
A' := XC\ (C, a). We have CCA'cA We claim that A = A'. 

Otherwise, as one can easily see, adding all the irreducibles that decom- 
pose A' to the family S, we obtain a nested family that contains properly 
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S. This contradicts the maximality of S. Clearly A = A' implies that 
(A) = (C) Co and thus dim(A) = dim(C) + 1. 

(3) We proceed by induction on s = dim(V). 

If s = 1 there is nothing to prove, there is a unique set complete and 
irreducible namely X. 

Let s > 1. Decompose X = Li^ =1 Xh into irreducibles. 

We have that a MNS in X is the union of MNS in each X;. Then s = 



Thus we can assume that X is irreducible. In this case we have that 
X € S for every MNS S. 

Let Bi, . . . , B s be the elements of S properly contained in X and maximal 
with this property. 

The set S consists of X and of the subsets Si := {A £ S \ A C Bi}. 

Clearly Si is a MNS relative to the set Bi (otherwise we could add an 
element to Si and to S contradicting the maximality of S). 

By induction Si has dim(Sj) elements and thus by (2) the claim follows. 



In the example of configuration spaces X spans the hyperplane in C s 
where the sum of coordinates equals to zero. We have described irreducibles 
via the corresponding subsets in {1, ... , s}. Under this correspondence, one 
can see easily that a maximal nested set S is formed by s — 1 elements and for 
any A € S with a > 2 elements either A contains a unique maximal element 
B € S with necessarily a — 1 elements or exactly two maximal elements 
B ± ,B 2 eS with A = B 1 U B 2 . 

We can present such a MNS in a convenient way as a planar binary rooted 
tree with s leaves labelled by {1, . . . , s}. Every internal vertex of the the 
corresponds to the set of numbers that appear on its leaves. 

For example the graph 



dim((X))=Et 1 dim((X,)). 



□ 



(62) 
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represents the MNS {1, 2}, {1, 2, 3}, {4, 5}, {1, 2, 3, 
Given a MNS S let us define a map: 



4,5}. 



PS ■ X -> S 
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as follows. Since Li Aes A = X every element a € X lies in at least an 
A € S. Also if an element a € X appears in two elements A, B € S the 
two elements must be necessarily comparable, thus there exists a minimum 
among the two. It follows that for any element a € X there is a minimum 
element ps{o) € S. 
Now a new definition: 

Definition 13.11. We shall say that a basis b := {a±, . . . ,a s } C X of C s , 
is adapted to the MNS S if the map a% i— » ps{ai) is a bijection. 

Such a basis always exists. It suffices to take, as in the proof of ll3.1Ul for 
every yl € S an element a € A — Uj-Bj, where the -Bj are the elements of S 
properly contained in A. 

Given any basis b := {a\, . . . ,a s } C X, we shall build a MNS <S& to which it 
is adapted, in the following way. Consider for any 1 < i < s the complete set 
Ai := X n ({di, ...,a s }) = { ai ,...,a a }. Clearly A x = X D A 2 D ■ ■ ■ D A s . 

For each i consider all the irreducibles in the decomposition of Ai . Clearly 
for different i we can obtain also several times the same irreducible, in any 
case we have: 

Theorem 13.12. The family Sb of all the (distinct) irreducibles that appear 
in the decompositions of the sets Ai form a MNS to which the basis b is 
adapted. 

Proof. By induction. Decompose X = A\ = B\ U B 2 U • • • U into irre- 
ducibles, by construction: 

k 

s = dim(Ai) = y^dim(-Bj). 

i=i 

We have that A 2 = (A 2 nB 1 )U(A 2 nB2)U - ■ -U(A 2 nB k ) is a decomposition 
of A 2 , not necessarily into irreducibles. 
Since dim(A 2 ) = s — 1 we have: 

k 

s - 1 = dim(A 2 ) = ^dim(A 2 n 
i=l 

Therefore dim(^2 H Sj) < dim{Bi) for exactly one index iq. In other 
words we must have that A 2 f]Bi = Bi for all the i ^ i$ . For such an index 
necessarily a\ G Bi Q . 

By induction, the family of all the (distinct) irreducibles that appear in 
the decompositions of the sets Aj, i > 2 form a MNS for (A 2 ), with adapted 
basis {a 2 , . . . , a n }. To this set we must thus only add Bi in order to obtain 
Sb- Thus Sb is a nested set with s elements, hence maximal and the basis b 
is adapted. □ 

Remark 13.13. One can easily verify that, conversely, every MNS S is of the 
form Sb for some adapted basis. 
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13.14. Non linear coordinates. Now we pass to the fundamental geomet- 
ric construction. 

Given a MNS S and a basis b := {di, . . . , a n } adapted to S let us consider 
the system of linear coordinates on U. 

If ps(ai) = A we denote also Oj := a a- We build now new coordinates 
ZAi A G S using the monomial expressions: 

(63) a a ■= J} z B - 

bg5, acb 

Given 4 G S, let 5a := {-B C A, Be S}. Clearly S A is a MNS for A (in 
place of X)and the elements ag with B G 5^ form a basis of A, adapted to 
S A - 

Take a G X with ps(a) = A and write a = Y1b<=s a CbCLb > ° b ^ ^et us 
substitute now the expressions (jfi3*|) : 

o-a ■■= ^ ° b n z ° = 

Bes A ces, bcc 

(64) H z B {c A + Yl Cb II 

Be5, AC-B BeS A , B ¥"A C&S, BCCCA 

Since A is the minimum set of S containing a we must have ca ^ 0. 

At this point we can proceed as in Section fl 1 . 51 and define an embedding 
of Rx into the ring of Laurent series in the variables za, A G S. Thus for a 
top form ip we can define a local residue ressfiip ( or ress4> if k is clear from 
the context). 

13.15. Proper nested sets. We are going to tie the concepts of the pre- 
vious section to that of no broken bases. Our goal is to define a bijection 
between NB(X) and a suitable family of MNS which we shall call proper 
nested sets. We need some preliminary steps. 

Assume thus that the basis b := {ai, . . . , a n } C X is no broken, we get: 

Lemma 13.16. a\ is the minimum element of pg b (ai) for every i. 

Proof. Let A := ps b (ai), by the definition of Sb, we must have that A de- 
composes one of the sets A^ = (a^, ...,a n ) D X. Necessarily it must be 
k <i. 

On the other hand belongs to one of the irreducibles of A{ , that there- 
fore is contained in each irreducible B of Ak, k < i that contains a,. 

It follows that A must be one of the irreducibles decomposing Ai. By defi- 
nition of no broken basis, a, is the minimum element of Ai = (oj, . . . , a n ) fl X 
hence also the minimum element of A. □ 

This property suggests us to define: 

Definition 13.17. A MNS S is said to be proper, if the elements 
as := mina, a G S \ S G S, form a basis. 
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Lemma 13.18. If S is proper, ps(as) = S. Thus the elements as ■= 
min a, a € S \ S £ S form a basis adapted to S. 

Proof. Let U € S with as £ U . If U C S, then as = min a G £/ and thus 
as = at/- 



If S is proper, we order its subsets Si, . . . ,S n using the increasing order 
of the elements as- We then set a% := as v and we have: 

Theorem 13.19. (1) The basis b := {a\, . . . , a n } is no broken. 

(2) In this way we establish a 1-1 correspondence between no broken bases 
and proper MNS. 

Proof. (1) From the Remark 113.91 setting Aj := Uj>iSj, we have that Aj is 
complete and decomposed by the Sj which are maximal. 

Clearly, by definition, a% is the minimum of Aj. It suffices to prove that 
Aj = (aj, . . . , a s ) n X. Hence that (Aj) = (dj, . . . , a s ) since Aj is complete. 

We prove it by induction. If i = 1 the maximality of S implies that 
A\ = X so di is the minimum element in X. Now ai ^ A2 so A2 7^ X and, 
since A2 is complete dim((A2}) < s. 

Clearly (A 2 ) D (a 2 , . . . ,a s ) so (A 2 ) = (a 2 , . . . ,a s ). 

At this point it follows that, by induction, (02, . . . , a s ) is a no broken basis 
in A2 and thus, since a\ is minimum in X, (a\, 02, • ■ ■ , a s ) is a no broken 
basis. 

(2) From the proof, it follows that the two constructions, of the MNS 
associated to a no broken basis and of the no broken basis associated to a 
proper MNS, are inverse of each other and thus the 1-1 correspondence is 
established. □ 



We can now complete our analysis computing the residues. Recall that, 
if b := {bi, . . . , b n } C X is a basis, we denote with: 

uj k := dlog(hi) A • • • A dlog{b n ) 
Theorem 14.1. Given two no broken bases b,c we have: 



Proof. We prove first that ress b Ub = 1- 

By definition bi = H A . cB z B hence dlog(6j) = J^AiCB dlog(z B )- 
When we expand the product we get a sum of products of type d log^^ ) A 
d\og{z B2 ) A • • • A dlog(zB n ) with Aj C Bi. 

Now, the unique non decreasing and injective map, of in itself is the 



Since the elements as are a basis, this implies that S = U. 



□ 



14. Residues and cycles 




identity. 
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Therefore in this sum, all the monomials vanish except for the monomial 
dlog(zi) A c£log(z2) A • • • A dlog(z n ), that has residue 1. 

Let us pass now to the second case b^c. This follows immediately from 
the following Lemma. □ 

Lemma 14.2. 1. If b ^ c, the basis c is not adapted to Sb- 

2. If a basis c is not adapted to S = {Si, . . . , £*„}, a MNS, we have 
resstOc = 0. 

Proof. (1) Let c = {c\, . . . , c„} be adapted to <S&, we want prove that we 
have b = c. We know that c\ = a\ = b±, is the minimum element of X. 

Let A be the irreducible component of X containing a\. This, by defini- 
tion, is an element of Sb- We have thus A = S±. We claim that ps b (ai) = A. 

This follows from the fact that a\ = b\ and Sb is proper. 

Set X' := S2 U • • • U S n , X' is complete. The set S2, ■ ■ ■ , S n coincides with 
the proper MNS S b ' associated to the no broken basis b[ := {62, • • • > b n } of 
(X'). 

Moreover clearly, d = {c2, . . . , c n } is adapted to S b '- Therefore b' = d by 
induction. Hence b = c. 

Using the first part, the proof of (2) follows the same lines as the proof 
of Lemma 111. 61 and we leave it to the reader. □ 

Remark 14.3. We end this section pointig out that by what we have proved 
it follows that for any no broken basis b we have that 

resb = res Sk 
with resb defined in Section fll.51 

The advantage of this new definition is that one sometimes uses mono- 
mial transformations of smaller degree and this could provide more efficient 
algorithms. 

14.4. A minimal model. Although we do not use it explicitly it may useful 
to understand the origin of the non linear coordinates we have been using 
and the entire theory of irreducibles and nested sets that we have built in 

Start from a family of hyperplanes in U = V*, given by a list X C V of 
linear equations, with complement Ax C U. 

In |19j we construct a minimal smooth variety Zx containing Ax as an 
open set with complement a normal crossings divisor, plus a proper map 
7r : Zx — » U extending the identity of Ax ■ 

The smooth irreducible components of the boundary of Ax are indexed by 
the irreducible subsets. To describe the intersection pattern between these 
divisors, in JH] we developed the general theory of nested sets. 

Maximal nested sets correspond to special points at infinity, intersections 
of these boundary divisors. In the papers [HI] and implicitly the authors 
use the points at infinity coming from complete flags which correspond, in 
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the philosophy of JH]) to a maximal model with normal crossings. It is 
thus not a surprise that by passing from a maximal to a minimal model the 
combinatorics gets simplified and the constructions become more canonical. 

Let us recall without proofs the main construction of [19j . 

For each irreducible S C X we have an orthogonal subspace S 1 - C U 
where S 1 - = {a G U | x(a) = 0, Vx G S}. 

From the collections of projective spaces P(J7/S ,_L ) of lines in U/S ± we 
deduce a map i : Aa — > U* x^gj P([//5" L ). Set Zx equal to the closure of 
the image i(Ax) in this product. In ^H] we have seen that Zx is a smooth 
variety containing a copy of Ax and the complement of Ax in Zx is a union 
of smooth irreducible divisors D$, having transversal intersection, indexed 
by the elements S G X. 

Still in |19| we showed that a family Ds i of divisors indexed by irreducibles 
Si has non empty intersection (which is then smooth irreducible) if and only 
if the family is nested. In particular a maximal nested set J\f identifies a 
special point at infinity pj\f, intersection of the s = dim({7) divisors corre- 
sponding to the irreducibles in J\f. The non linear coordinates are indeed 
coordinates in a local chart around pj\f in which pjj = and the boundary 
divisors are given by the vanishing of the s coordinates. 

15. Final considerations 

As we have seen, the actual problem of computing explicitly the functions 
which we have discussed can be approached through several different ways. 
We have proposed three approaches, by expansion into partial fractions, 
by solving a system of linear equations interpreting the defining differential 
equations and finally with the method of residues. It is not really clear to 
us which is the most efficient. 

A lot of computer work, aiming at computing Clebsch-Gordan coefficients 
using these methods has been done by several authors using the method of 
residues but whether this is really the fastest algorithm is yet unclear. 

Due to limited space we have discussed few examples. For our motivations 
one of the most interesting examples is the set X of positive roots of a given 
root system. There are still many things to be uncovered in this case. 
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